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Abstract 

Let h(x, y) be a non-degenerate binary cubic form with integral coeffi¬ 
cients, and let S be an arbitrary finite set of prime numbers. By a classical 
theorem of Mahler, there are only finitely many pairs of relatively prime 
integers x,y such that h(x,y) is an S-unit. In the present paper, we re¬ 
verse a well known argument, which seems to go back to Shafarevich, and 
use the modularity of elliptic curves over Q to give upper bounds for the 
number of solutions of such a Thue-Mahler equation. In addition, our 
methods gives an effective method for determining all solutions, and we 
use Cremona’s Elliptic Curve Database to give a wide range of numerical 
examples. 
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1 Introduction 

Let h(x,y) be a non-degenerate cubic form with integer coefficients, and let S 
be a finite set consisting of s distinct prime numbers, say pi,P 2 , ■ • • ,p s - Then 
the Thue-Mahler equation 


h{x,y) = ±Y[pT (1) 

i=i 

has finitely many solutions among relatively prime integers x , y and non-negative 
integers ei, e^, ■ ■ ■ , e s . In geometric terms, if we denote by Zg the ring of S- 
integers, and by Y the affine variety defined as the complement of zeros of h(x, y) 
in a projective line, then the solutions of above Thue-Mahler equations, modulo 
the identification of (x,y) and (—x,—y), bijectively correspond to the elements 
of Y(Zg). 

Mahler gave an ineffective proof of the finiteness of Y(Zg), and Coates [5], 
0 , m later obtained an effective finiteness of Y(Zg) using Baker’s estimate of 
linear forms in logarithms together with its p-adic analogues. However, explicit 
determination of Y(Zg) based on Baker’s method is often practically impossible 
due to the astronomical size of resulting upper bound for the height of a putative 
solution t £ Y(Zg). 

The aim of the current article is to present a new approach to Thue-Mahler 
equation. In order compute Y(Zg), we design a descent procedure, which mim¬ 
ics the Kummer homomorphism for rational points on elliptic curves. More 
precisely, we will construct a natural map 

k : Y{Zg) —> {Elliptic Curves over Q up to isomorphism} (2) 

t ► X t (3) 

which associates an elliptic curve X t to an unknown solution t £ Y(Zg), and 
study local properties of X t . In particular, we will show that X t has good 
reduction outside of S , discriminant of h(x 7 y ), and 2. It allows one to compute 
the image of k without knowing elements of Y(Zg). On the other hand, for 
an elliptic curve E, we will show that k~ x {E) is naturally a zero dimensional 
algebraic variety defined by explicit polynomials with rational coefficients, whose 
Q-points correspond t £ Y( Zg) equipped with an isomorphism from X t to 
E. In particular, one can numerically compute n^ 1 (E) from coefficients of a 
Weierstrass equation for E. 

Existence of such a map k first allows us to bound the cardinality of Y(Zg) 
from above, in terms of the number of elliptic curves whose conductor belongs to 
a finite list of integers, where the list of possible conductors is obtained from the 
coefficients of h(x , y) and S. The number of such elliptic curves can be bounded 
from above either using the work |T| of Brumer and Silverman or modularity of 
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elliptic curves. The former has better asymptotics, while the latter provides a 
practical algorithm. 

Theorem 1.0.1. Let S be a finite set of primes containing 2 and prime divisors 
of the discriminant ofh(x,y). Let G(S) be the number of isomorphism classes 
of elliptic curves which have good reduction outside of S. Then we have 

Y(Z S ) < |AutQ(F)| x G(S). (4) 

Combining it with an upper bound for G(S), due to Brumer and Silverman, we 
have 


Y(Z S ) < |Autnj(y)| x fc 2 M5 +£ . (5) 

where M is the product of all prime numbers in S, e is an arbitrary positive 
number, and k 2 is a constant depending on e. 

In fact, one can identify AutQ(T) with a subgroup of the symmetric group 
acting on zeros of h(x,y), so it has at most six elements. 

We would like to stress that our proof is manifestly constructive, which ulti¬ 
mately relies on the modularity of elliptic curves defined over rational numbers. 
More precisely, we give an explicit characterisation of fibres of k which allows 
us to compute Y(Zg) from k (Y(Z,g)), and the modularity of elliptic curves pro¬ 
vides a constructive finiteness of re(Y(Z,g)). In order to show that our approach 
to compute Y(Zg) works in practice, we append the tables of complete solutions 
for the following equations. 


h(x,y) 

S 

Table 

x(x - y)y 

{2,7,11,13} 

m 

x(x - y)y 

{2,3,431} 

m 

x(x - y)y 

{2,3,5,53} 


(x 2 + 7 y 2 )y 

{2,3,5, 7} 

si 

C x 2 + 7 y 2 )y 

{2,7,11,13} 

si 

C x 2 + 3 y 2 )y 

{2,3,11} 

El 

2{x 2 +y 2 )y 

{2,3,7,11} 

□ 

(x 2 + y 2 )y 

{2,5,13} 

M 

(■ x 2 - 2 y 2 )y 

{2,7,29} 

M 

{x 2 - 2 y 2 )y 

{2,7,29} 

m 

( x 2 - 3 y 2 )y 

{2,5,7,11} 

sn 

{x 2 - 7y 2 )y 

{2,3,7,11} 

\m 

x 3 — x 2 y — 4 xy 2 — y 3 

{2,5,13} 

S3] 

x 6 — x 2 y — 2 xy 2 — 2 y 6 

{2,5,19} 

si 

x 3 + y a 

{2,3,5} 

m 

x J + 2y i 

{2,3,5} 

ss 

x a - y 3 

{2,3,5} 

S3 

x : 3 - 2y i 

{2,3,5} 

S3 
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Note that we omitted the trivial solution h( 1,0) = 1 in the appended tables, 
and Table [3] had to be abbreviated due to a large number of solutions. 

The choices of h(x,y ) and S' in © had to be restricted according to our 
computational capability, and the choices are made to show the flexibility that 
we have. The main restriction comes from one’s ability to find the C 4 and 
Cq invariants of all elliptic curves whose conductor divides the worst possible 
conductor given in Proposition 12.3.11 Other computational difficulties are neg¬ 
ligible. Computational issues are further discussed in Section [7] 

When we computed the solutions of the equations listed above, we exploited 
Cremona’s Elliptic Curve Database, from which we read the coefficients of el¬ 
liptic curves with specified conductor. After that, we compute k~ 1 (E) for each 
curve E read from the database. Of course, it is highly non-trivial task to 
establish a complete list of isomorphism classes of elliptic curves of specified 
conductor, and we are outsourcing this job to Cremona. We note that this job 
is computational infeasible without modularity, and even with modularity it 
takes significant further efforts to obtain practically efficient algorithm. Never¬ 
theless, the absence of modularity is the main theoretical and practical obstacle 
to generalising our method to number fields. Once we have the coefficients of 
necessary elliptic curves, then it takes no more then 20 seconds to generate each 
of the tables listed above. 

The spectacular resolution of Fermat’s Last Theorem by means of modular 
methods as well as ID. 0, 0 and [5] uses level lowering argument in a crucial 
way, which allows one to produce obstruction to existence of a solution by 
numerically showing that certain space of modular forms is zero dimensional. 
In contrast, we will be using modularity theorem in order to produce a complete 
set of solutions for a given Thue-Mahler equation, without a priori guess on the 
number of solutions. In particular, we do not use any form of level lowering. 

We outline the contents of the paper. In Section [2j we define k and study 
basic properties of the elliptic curve X t associated to £ G Y(Zs). In Section [3l 
we study n~ l (E) for an elliptic curve E given in terms of a Weierstrass equa¬ 
tion. In Section [TJ we prove the main theorem on the upper bound of the 
cardinality of Y(Zs), compare our upper bound with Evertse’s upper bound. 
The proof is manifestly constructive, and it provides an algorithm to determine 
Y(Zs). We discuss the algorithmic aspect in Section [5] In Section [ 6 l we ex¬ 
plain how the algorithm is implemented in the computer algebra package, and 
discuss its performance. We also list the cardinalities of Y(Zs) as we vary S. 
In Section [71 we recall the work of Tzanakis and de Weger which proposed, 
based on Baker’s method and further optimisations, a practical algorithm for 
Thue-Mahler equations, and we compare the natures of two approaches. In 
Section[Sl we specialise h(x, y) in order to explain the connection to generalised 
Ramanujan-Nagell equations. Some particular generalised Ramanujan-Nagell 
equations are solved, from which we observe a pattern among the number of 
solutions as we vary S. 
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2 Definition of k and its properties 

Consider a binary cubic form 

h(x, y ) = ax 3 + bx 2 y + cxy 2 + dy 3 (7) 

with relatively prime integer coefficients. We assume that the discriminant 

S = 3b 2 c 2 + 6abcd — 4 b 3 d — 4ac 3 — a 2 d 2 (8) 

of h{x, y) is non-zero, or equivalently that h(x , y) has three projectively non¬ 
equivalent zeros over an extension of Q. 

Let P* be the projective line with homogeneous coordinates x and y. Let Z 
be the subscheme of P \ y defined by h(x, y) = 0, and let Y to be the complement 

Y = P 1 - Z (9) 

which we view as an affine variety embedded in P), y . In particular, a point t in 
Y(R) for some ring R will be represented as a pair (x t : yt) £ P 1 (i?) such that 
h{xt,yt ) is a unit in R. 

The aim of this section is to introduce the map k which associates an ellip¬ 
tic curve to a point in Y(R), and to study its basic properties. We will first 
construct a generically smooth map X —> Y, and the associated elliptic curve 
will be constructed by taking its fibre. The associated elliptic curve is naturally 
equipped with additional structures, which we will analyse in this section. 

2.1 Coordinates of Y 

Our definition of Y as an open subscheme of P 3 y , endows Y with homogeneous 
coordinates x and y. but we would like to introduce another coordinate e of Y , 
which makes our later discussion simpler. Consider Y C A 3 ye , defined by 

Y:h(x,y)e = 1 (10) 

where h(x, y) is the defining equation of Z = P 1 — Y. Let G m act on A 3 ye with 
weight 1, 1, and —3. That is to say, for any ring R, A € R x , and (x,y,e) £ 
A xye(R)y the action of A is given by 

A • (x,y,e) = (Ax, Ay, X~ 3 e) . (11) 

Because the action preserves GOD, one can consider the quotient G m \Y, which 
is just Y. Indeed, x and y are homogeneous coordinates of degree one, defining 
the embedding G m \Y —> P 1 . The coordinate e of Y is redundant, but it will be 
convenient for later purposes. 
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2.2 Construction of the family /: X —> Y. 

We describe /: X —> Y in this subsection. We will define X as a quotient of X, 
where X is an affine subscheme of A 3 ye x A^ vw . The defining equations of X 
are 


h(x,y)e = 1 (12) 

e ■ w 2 = h(u, v)(yu — xv). (13) 

Now we let G m x G m act on X in the following way. If R is a ring, (A, /i) G 
G m {R) x G m (R), and (x,y,e,u,v,w) G X, then we define 

(A, y) • (x, y, e, u, v, w ) = (Ax, Ay, A~ 3 e, yu, yv, X 2 y 2 w). (14) 

Since m and m are preserved by the action m, we may define 

X = G m x G m \X. (15) 

Furthermore, the projection X —> Y descends to X —> Y, which we denote by 

/■ 

Remark 2.2.1. Geometrically speaking, X parametrises the double covers of 
P 1 branched along the divisor Z U {f} of degree four, as t varies in Y. However, 
this does not uniquely characterise X , since there are more than one such double 
covers which are not isomorphic to each other over Q. 

Remark 2.2.2. When h(x, y) = x(x — y)y, Y can be identified with the affine 
line without 0 and 1, by taking A = x/y as an affine coordinate. Then, the 
equation 

w 2 = u(u — v)v(u — Xv) (16) 

defines the Legendre family of elliptic curves over Y. Our family X —> Y of 
elliptic curves in this case is represented by 

(x(x — y)y 2 )~ lr w 2 = u(u — v)v(u — Xv) (17) 

so one can view it as a quadratic twist of the original Legendre family by x(x — 
y)y 2 . This relation between Legendre family and our X is available because of 
the affine coordinate A for Y. When h(x, y) has no rational linear factor, such 
an affine coordinate is not available, whence the original Legendre family does 
not directly generalise. It is the advantage of the twisted family (ED that it 
generalises to general h(x,y). 

2.3 Properties of /: X —> Y 

In this subsection, we study basic properties of /: X Y. For some ring R 
and t = (t x : t y ) G Y(R), the elliptic curve X t is defined by the equation 

X t : e • w 2 = h(u, v)(y t u - x t v) (18) 
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which be interpreted as a quadratic twist of 

X' t : w 2 = h(u,v)(y t u - x t v) (19) 

by e, although X' t does not patch together to form a family of elliptic curves 
over Y. In fact, X[ is not even well-defined on the projective equivalence class 
of t = (xt : yt), and only its quadratic twist X t is well-defined. In any case, we 
can compute the discriminant of the right hand side of (1191) , as well as prove its 
properties. 

Proposition 2.3.1. The discriminant of right hand side of (1191) is 

h{x t ,yt) 2 ■ 5 (20) 

where 5 is the discriminant of h{x,y). Furthermore, we have: 

1. if t £ Y(Zs), then X' t and has good reduction outside of S and 25. 

2. if an odd prime p £ S is a prime of bad reduction for X' t and p does not 
divide 25, then h{x , y) = 0 has at least one solution modulo p. 

3. if p does not divide 25, then X' t has either good or multiplicative reduction. 

Proof. The formula (1201) follows from representation of the discriminant in terms 
of differences of roots. Indeed, if P(x) is a polynomial of one variable with roots 
,a n , then the discriminant 5p of P{x) is 

J^( a * — a j ) ■ (21) 

i<j 


If P{x) = (x — fl)Q(x), and a n = /?, then the above formula can be rewritten 

P 08 -«<) a x II ( 22 ) 

i<n i<j<n 


which equals Pi(/?) 2 ■ 5q. 

If t £ Y(Zs), and (x t : yt) is some representative of t, then X[ has good 
reduction away from S and 25. Indeed, a double cover of P 1 branched along 
four distinct points is smooth away from characteristic two. 

Suppose p £ S does not divide 2<5 and X' t has bad reduction at p. Since 
p does not divide 5 , h(x,y) has three distinct roots modulo p. Thus X' t has 
bad reduction if and only if t coincides with one of three zeroes of h(x,y). In 
particular, t is a solution of h(x, y) = 0 modulo p. 

Assume that p does not divide 25 and X) has a bad reduction at p. Then the 
right hand side of (fl9l) cannot have a cubic factor, since such a factor will force 
h(x, y) to have at least a square factor modulo p, contradicting the assumption. 
In other words, X' t has either good or multiplicative reduction. This completes 
the proof of the proposition. □ 

For Xt, we can prove the following:- 
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Proposition 2.3.2. With the notation as in the previous proposition, the dis¬ 
criminant of m is 

h{x t ,y t ) 4 5. (23) 

In particular, for any t £ Y(fLg), X t has good reduction outside of S and 25. 

Proof. If follows immediately from Proposition 12.3.11 and the description of the 
discriminant in terms of differences of roots. □ 

Remark 2.3.1. The dependence of X[ on the choice of representative (xt ■ yt) 
is not so serious, as far as we work with rational numbers. We can always take 
(xt ■ yt), such that Xt and y t are relatively prime integers and x t is non-negative. 
If we work over a number field which has either positive class number or more 
units than ±1, this is not straightforward. Working with X' t has advantage that 
the conductor of X' t is usually smaller than X t and it has in some sense finer 
information about t than X t does. On the other hand, X t is associated naturally 
from t, whose isomorphism class is independent of choice of representative for 
t, so it is technically more convenient. 

3 Defining equations of fibres of k 

In the previous section, we defined a map 

k '. 1 1 —y X t (24) 

which associates an elliptic curve X t to a solution t £ Y(Zs). The aim of the 
present section is to describe k~ 1 {E) as a variety defined by explicit polynomials. 

3.1 Some invariant theory 

We briefly review basic invariant theory of binary forms which are relevant for 
us. We start with invariants of binary quartic forms. 

Let 

q = A 0 u 4 + Aiuv 3 + A 2 U 2 v 2 + A 3 uv 3 + A^v 4 (25) 

be a generic binary quartic form in u,v, with coefficients Aj’s. We choose 

h = —^2 ~~ ~^AiA 3 + ^ 0^4 (26) 

h = 2l6 Ai - + T6 A ° A * + T6 A ' A4 - 

as generators for the ring of invariants of binary quartic forms. Note that they 
have rational coefficients, and their degrees are two and three respectively. 

In fact, these two invariants are algebraically independent, or the ring of 
invariants is isomorphic to the polynomial ring in two variables. To stress the 
dependence of the invariants on the coefficients, we denote by q(A) the quartic 
form with coefficients A = (Ao, A\, A 2 , A 3 , A±), and denote their invariants by 
hi A) and h{A) respectively. 
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Proposition 3.1.1. Let q(A) and q(A') are two binary cubic forms with rational 
coefficients. They are linearly equivalent over Q if and only there exists A £ Q x 
such that 

h(A) = A 2 (A') (28) 

h(A) = A 3 04') (29) 

holds. 

Proof. Classical invariant theory. □ 

Now we consider the invariant theory of a pair of binary forms. Let 

q(A , B) = (Aqu + Aiv)(Bqu 3 + Biu 2 v + B 2 UV 2 + B 3 v 3 ) (30) 

be a product of a binary linear form and a binary cubic form, where A = (Alo, A\) 
and B = (Bo, • • • , Bo) denote the coefficients of a linear form and a cubic form 
respectively. Two invariants 

c 4 (A, B)=- 16(- A\b\ + 2,AlBoB 2 + AoA iBiB 2 - A 2 B‘ 2 - 9A 0 A 1 B 0 B 3 + 3 A^BiB z ) 

(31) 

c 6 (A,B) = - 32(2 A\Bl - 9A(B 0 B l B 2 - 3 A 0 A 2 B^B 2 + 18 A 0 A 2 B 0 B^ - SA^B^ 

(32) 

+ 2A 3 B 3 + 27 A 3 B 2 B 3 - 27A 0 A 2 B 0 B 1 B 3 + - 27 A 2 A 1 B 0 B 2 B 3 

(33) 

- 9A 3 B 1 B 2 B 3 + 27A 3 B 0 B 2 3 ) (34) 

generate the ring of invariants. 

We digress for a discussion on the relation between above invariants and 
invariants of an elliptic curve often used in the literature. If an elliptic curve E 
is given by 

E : y 2 = x 3 + a 2 x 2 + a 4 x + a 6 (35) 

then 

c 4 (E) =1602 — 48a4 (36) 

ce(E) = — 64a2 + 2880204 — 864o6 (37) 

are often called C4-invariant and C6-invariant of E. If we take Aq = 0, A\ = 1, 

Bq = 1, Bi = 02, B 2 = 04, and B 4 = ae, then c 4 (A,B ) and cq(A,B ) are 
precisely the C 4 -invariant and C6-invariant of the elliptic curve defined by (1551) . 

The discriminant of (1351) is given by 

6(E) = —16(— a^a 2 + 4a|a6 + 4a| — 18020406 + 27og) (38) 

and it can be alternatively written as 

6(E) = 2 6 3 3 (c 4 (E) 3 - c 6 (E) 2 ) = 1728 (c 4 (E) 3 - c 6 (E) 2 ). (39) 
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The discriminant of q{A,B) viewed as a single binary quartic form, its discrim¬ 
inant is 

S(q(A, B)) = - (- B\B\ + AB 0 B 3 + 4B?B 3 - 1 ^B n B x B 2 B z + 27 B 2 B 2 ) (40) 
x (— A^Bq + AqA^Bi — AqAiB2 + AqBz ) 2 . (41) 

If we take A$ = 0, Ai = 1, Bo = 1, B 4 = 02 , -B 2 = 04 , and B 4 = ae, then two 
discriminants are related by 


16 ■ 6 (q(A,B)) = 6 (E). (42) 

The relation between C4, cq and the previously introduced I 2 and I 3 are more 
straightforward. Indeed, they are related by 

d(q(A,B)) =192- I 2 (q(A, B)) (43) 

=2 6 3 • I 2 {q{A, B)) (44) 

c 6 (q(A, B)) = - 13824 • I 3 (q(A, B)) (45) 

= -2 9 3 3 -I 3 (q(A 1 B)) (46) 


where we view q(A , B) as a product of two forms on the left hand side, while 
on the right hand side we view it as a single quartic form whose coefficients are 
quadratic forms in Aj's and Bi s. 

Now we return to the invariant theory of q(A,B). 

Proposition 3.1.2. Let q(A , B) and q(A', B') are two binary quartic forms with 
factorisation as a product of a linear and a cubic factor, and suppose that the 
coefficients A,A',B and B' are rational numbers. They are linearly equivalent 
over the rational numbers if and only if there exists A £ Q x such that 

c 4 (q(A,B))=\ 4 C 4 (q(A',B')) (47) 

c 6 (q(A,B))=X e c 6 (q(A / ,B 1 )) (48) 

holds. 

Proof. Classical invariant theory. □ 

3.2 Faithfulness of descent 

The overall strategy is to study t £ Y(Zg) in terms of X t . In other words, we 
consider the map 

n \ 1 1 —y Xt (49) 

and try to use k in order to compute Y(Zg). To realise this, we will show two 
key properties of k: 

1. K is n-to-1 map, where n is an explicit integer less than six. 
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2. Given a Weierstrass equation of an elliptic curve E, one can compute 
n~\E). 

Let us consider the first property of k. Let E be an elliptic curve. We would 
like to count the number of t £ Y(fLs) for which Xt is isomorphic to E, where 
E is given as a Weierstrass equation 



E : y 2 + a\y + 03 = x 3 + a- 2 X 2 + a^x + ae 

(50) 

with rational coefficients. On the other hand, X t is defined by 



ew 2 = h{u , v)(yu — xv ) 

(51) 

where h(x , y)e 

= 1. We rewrite X t as 



w 2 = h(u, v)(yu — xv)h{x, y) 

(52) 

and let 

Q{x, y, h) = (yu - xv) ■ h(u, v)h(x, y) 

(53) 


be the associated joint quartic form in u and v. In fact, one could write 
Q(x,y,h) = Q(t,h), in the sense that Q(Xx,Xy,h) is rationally equivalent to 
Q(x,y,h) for any A £ Q x . For simplicity of notation, let 

Ci{x, y, h) = c 4 (Q( x, y, h)) (54) 

c 6 (x,y,h) = c 6 (Q(x,y,h)) (55) 

be the invariants of Q(x, y, h). The set of all t £ Y(Zs) for which X t is ismorphic 
to E is defined by the equations 

C 4 ,{x, y, h) = X 4 ca(E) (56) 

c e (x, y, h) = X 6 c e (E) (57) 

where a(x,y,h) and ce(x,y,h) are homogeneous polynomials of degree eight 
and twelve in x,y, respectively. In fact, we can eliminate A from above two 
equations to obtain 

J 24 :{x, y, h, E) := A 12 (c 6 (E) 2 C 4 {x, y, h ) 3 - ca(E) 3 c & (x, y, h) 2 ) (58) 

which factors as 

J- 24 { X, y, h, E ) =A 12 • h(x, y ) 6 ■ J 6 ( x, y, h, E) (59) 

where Jq{x , y 1 h, E) is the homogeneous polynomial of degree six in variables x 
and y, characterised by above equality. 

Proposition 3.2.1. The set of points t £ Y(Zs) for which X t is isomor¬ 
phic to E is in bijection with the projective equivalence class of solutions of 
Je(x,y, h, E) = 0 such that (l56l) and (l57l) have a common solution in X. In 
particular, this set has cardinality at most six. 
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Proof. We first show that Je(x,y,h,E), viewed as a homogeneous polynomial 
in x and y, is not identically zero for an elliptic curve E and a non-degenerate 
binary cubic form h(x,y). We work with complex numbers, although any alge¬ 
braically closed field of characteristic zero suffices our purpose. Let (xj, yt) with 
i = 1,2, • • • , r be a sequence of non-equivalent complex zeros of Je(x,y,h,E). 
They are precisely the values for which 

Qi(xi,yi,h) = (ym - Xiv) ■ h{u,v)h{xi,yi) (60) 

becomes equivalent to 

Qe '■= v(u 3 + ci 4 (E)uv 2 + (Iq(E)v 3 ) (61) 

as joint binary quartic forms in u and v. Isomorphisms between two joint forms 
Qi and Qe correspond to the projective isomorphisms between two divisors 
represented as zeros of Qi and Qe, which sends {xi,yf) to (0,1). Since E is 
an elliptic curve, divisor of u 3 + 04 (E)uv 2 + ae(E)v 3 consists of three distinct 
points. There are six projective automorphisms of P 3 y which sends the divisor 
of h(u, v) to the divisor of u 3 + cii{E)uv 2 + ae(E)v 3 . Thus r is at most six, and 
Je(x,y, h, E) cannot be identically zero. 

If t = {x : y) G Y (Zg), then h(x, y) ^ 0. Thus the vanishing of J 2 a{x, y, h, E) 
is equivalent to vanishing of Je(x,y,h, E), from factorisation of J 24 ,(x, y, h, E). 
If t further satisfies the condition X t being isomorphic to E , then clearly (1561) 
and m have a common solution. Conversely, if Je(x, y, h, E) vanishes at some 
point t, then for each solution satisfying (l56ll and d5ZD- Thus such t together 
with A gives rise to X t equipped with an isomorphism to E. This completes the 
proof. 

□ 

Corollary 3.2.1. For a fixed elliptic curve E, the number of t £ Y(Zg) for 
which X t is isomorphic to E is at most the cardinality o/AutiQ>(Y). 

Proof. It is implicit in the proof of Pronosition l3.2.n observing that the isomor¬ 
phisms between h(u, v ) and it 3 +ai{E)uv 2 +ae(E)v 3 is a torsor for Autjj(F). □ 

Before we move on, we analyse coefficients of Je(x,y, h, E). As we men¬ 
tioned earlier, it is homogeneous of degree six in x and y. With respect to the 
coefficients of h(x,y), namely a,b,c, and d, it is homogeneous of degree six as 
well. In terms of coefficients of E , it has degree twelve in the following sense. If 
we write E as 


y 2 + a\xy + a^y = x 3 + a^x 2 + 04 a; + ae (62) 

then J 6 (x,y, h, E) is a polynomial in variables 04 , 02 , 03,04 and Oq. If we take 
degree of a m to be to, then Jg(x,y,h,E) is homogeneous of degree twelve in 
o m ’s. More concretely, if we take a model of E for which 04 = 02 = 03 = 0, 
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then each term of Je(x,y, h, E), viewed as a polynomial in a 4 and as, is either 
a 3 or (Lq. Based on this observation, we arrange 

6 

Mx,y,h,E) = ]T (' C*{h)a\ + D*(h)a 2 ) x^Y (63) 

i=0 

where C* (h ) and D (/() are homogeneous polynomials of degree six in a, b, c, 
and d. In fact, all of them have a large integer factor, so we let 

C*(h) = 2 22 3 3 Ci(h) and D*(h) = 2 22 3 3 AW (64) 

and give formulas for Ci(h) and Di(h): 

C 0 (h) = (2c 3 - 9 bed + 27 ad 2 ) 2 
D 0 (h) = —27(—c 2 + Sbd ) 3 

C\(h) = —6(2 c 3 — 9 bed + 27 ad 2 )(—bc 2 + 6 b 2 d — 9 acd) 

Di(h) = — 81(— bc+ 9 ad)(—c 2 + 3 bd ) 2 

C 2 {h) = -3( 6 2 c 4 - 24ac 5 + 18 b 3 c 2 d + 90 abc 3 d - 108 b 4 d 2 + 216 ab 2 cd 2 - 567 a 2 c 2 d 2 + 486 a 2 bd 3 ) 
D 2 (h) = —81(—c 2 + 3bd)(2b 2 c 2 — Sac 3 — 3 b 3 d — 9abcd + 81 a 2 d 2 ) 

C 3 {h) = -2(13 b 3 c 3 - 72 abc 4 - 72b 4 cd + 567 ab 2 c 2 d - 432a 2 c 3 d - 432 ab 3 d 2 + 2ASa 2 bcd 2 + 729 a 3 d 3 ) 
D 3 (h) = —27 (—bc+ 9ad)(7b 2 c 2 — 18ac 3 — 186 3 d + 3 Qabcd + 81a 2 d 2 ) 

Ci{h) = —3 ( 6 4 c 2 + 18a6 2 c 3 — 108a 2 c 4 — 24& 5 d + 90 ab 3 cd + 216 a 2 bc 2 d — 567 a 2 b 2 d 2 + 486 a 3 cd 2 ) 
Di(h) = — 81(— b 2 + 3ac)(26 2 c 2 — 3ac 3 — 3 b 3 d — 9 abed + 81 a 2 d 2 ) 

05 (h) = 6(b 2 c — 6ac 2 + 9abd)(2b 3 — 9 abc + 27a 2 d) 

05 (h) = — 81(— bc+ 9 ad)(—b 2 + 3ac) 2 
C 6 (h) = (2b 3 - 9abc+27a 2 d) 2 
D 6 (h) = —27(—b 2 + Sac) 3 

Above coefficients seem to admit certain invariant theoretic meanings, but for 
the purpose of the current article, we only need to regard Je(x,y,h,E) as a 
defining equation of n~ l (E) given as explicitly as possible. 


4 Main theorems 

Now we are ready to state and prove main theorems. Before we state our result, 
recall that the strategy was to consider 

k: t*-+Xt (65) 

and study t in terms of X t . In the previous subsection, we focused on the 
analysis of k~ 3 (E) for a fixed E. What remains is to consider the image of k. 

Among several approaches to analyse the image of k, the common funda¬ 
mental property of k is that X t has good reduction outside from 26 and S. It 
is convenient to introduce a notation for the number of such curves. 
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Definition 4.0.1. Let N be a positive integer. Let <7 (IV) be the number of 
isomorphism classes of elliptic curves whose conductor is N. For a finite set S 
of primes, let 

M = p (66) 

pgs 

and let G(S) be the number of isomorphism classes of elliptic curves which have 
good reduction outside of S. 

Theorem 4.0.1. Let e be any positive integer. There exists constants k\, k 2 
depending on e such that 

g(N) < hN^ +£ (67) 

G(S) < k 2 M? +e ( 68 ) 

holds. If we assume (a part of) BSD conjecture and generalised Riemann hy¬ 
pothesis for the elliptic curves y 2 = x 3 + n, n £ Z, then there exists constants 
ki and k$ for which 

G(S) < kiM'°s iog*f (69) 

Proof. It is a theorem of Brumer and Silverman. □ 

As a corollary of above theorem, we obtain a first, rather crude, upper bound 
for the cardinality of Y(Zs). 

Theorem 4.0.2. A ssume that S contains 2 and prime divisors of S. Then 

\Y(Z S )\ < |Aut Q (Y)| x G(S). (70) 

In particular, 

\Y(Z S )\ < |Aut Q (Y)| x fe 2 M3 +e . (71) 

where M is the product of all primes numbers in S, and k 2 is a constant which 
only depends on an arbitrary positive number in Theorem \4-0.1\ Conditional on 
generalised Riemman hypothesis and BSD conjecture for the curve y 2 = x 3 + n, 
we have 

\Y(Zs)\ < |Autnj(Y)| x kiM^nku. ( 72 ) 

where k± and k$ are absolute constants. 

Proof. The first assertion follows from Proposition 13.2.11 and Corollary 13.2.11 
The second follows from the first assertion combined with Theorem 14.0. II □ 

Remark 4.0.1. In the course of deriving the above upper bound, we have 
forgotten all information of X t except the divisors of the conductors of X t . 
Numerical computations indicate that the reduction type at the primes dividing 
2 S takes a particular form, which will facilitate practical computations. 
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4.1 Comparison with Evertse’s bound 

Evertse proved the following remarkable upper bound:- 

Theorem 4.1.1. Let h(x,y) be an integral binary form of degree n > 3 which 
is divisible by at least three pairwise linearly independent linear forms in some 
algebraic number field. Let p\,p 2 , - • • ,Pt be a sequence of distinct primes. The 


equation 

t 

M*,v)\ = Hp? ( 73 ) 

i—1 

in relatively prime integers x,y and non-negative integers ki,k 2 , ••• ,kt has at 
most 

2 x 7 " 3 ( 2t + 3 ) ( 74 ) 

solutions. There is analogous explicit upper bound for number fields. 

Proof. Corollary 2 of [TO] □ 

As a direct consequence, we obtain 

|y(Z s )| < 7 9 < s+3) (75) 

when s is the cardinality of S. Let us take S to be the set 

S = {p : p < T} (76) 


of all primes up to a positive number T. Then by prime number theorem s 
is asymptotically T / log T. Therefore, Evertse’s bound can be rewritten, in a 
logarithmic scale, as 

log |r(Zs)| = o ( ■ 

On the other hand, the standard estimate shows that 

log M = ^2 lo §P = °( T ) 

p<T 

and our unconditional upper bound of Theorem 14.0.11 becomes 

log |E(Zg)| = O (T). 

The conditional upper bound of Theorem ll. 0.11 becomes 

which is comparable to Evertse’s. 

Putting aside the comparisons between upper bounds for the cardinalities 
Y(Z S ), we point out crucial difference of our method from Evertse’s in terms 
of effectiveness. Evertse’s upper bound is ineffective, in the sense that it does 
not provide an algorithm to decide T(Zg). In contrast, our proof is manifestly 
constructive, especially if one combines with modularity of elliptic curves. We 
elaborate on the constructive aspects of our method in the next section. 


(77) 

(78) 

(79) 

(80) 


15 









5 Algorithmic aspects: effective Mordell conjec¬ 
ture 

In this subsection, we elaborate on the algorithmic aspects of our proof, which 
answers the effective Mordell conjecture for Y. Effective finiteness of Y (Z 5 ) can 
be formulated in at least two ways: 

1 . to have an explicit upper bound on height of t € Y(fZg) in terms of 
coefficients of h(x, y) and S. 

2 . to have a procedure which enables one to determine Y(Zs), in provably 
finite amount of time, for a numerically given h(x, y) and S. 

The first version of effectiveness implies the second. Indeed, if T is obtained 
bound, then factoring h(m, n ) as m and n varies among all integers with ab¬ 
solute value at most T, one can determine Y(Zg). In fact, Baker’s method in 
principle provides such an upper bound. However, the efficiency of such proce¬ 
dure depends on the size of T, and the astronomical size of numbers obtained 
from Baker’s bound often prevents one from computing Y (Z 5 ) in practice. Our 
method directly provides the second version of effectiveness without a priori 
upper bound for height of t G Y(Zg), and we shall describe the procedure as 
it is implemented in a computer algebra package in order to generate tables of 
numerical examples. 

Recall that the principal tool for us is the map 

tv : t 1 — y Xf ( 81 ) 

which associates an elliptic curve X t to a putative solution t. Computation of 
k~ 1 (E) amounts to solving polynomials in one variable, such as Je(x,y,h,E). 
In particular, k~ 1 (E) can be effectively decided once the coefficients of E are 
known. Thus, what remains is to determine all possible elliptic curves E for 
which k~ 1 {E) is possibly non-empty. 

Theorem 5.0.2. For a finite set S of prime numbers, the coefficients of iso¬ 
morphism classes of elliptic curves which have good reduction outside of S can 
be determined algorithmically. 

Proof. We give a brief description. For details, especially practical issues, we 
refer to Cremona. Using modularity theorem, one may compute the isogeny 
classes of elliptic curves by means of modular forms. Space of modular forms 
can be computed using modular symbols for example. For each isogeny class 
of elliptic curves, one can decide the isomorphism classes of elliptic curves con¬ 
tained in it. In fact, there are at most eight isomorphism classes in a fixed 
isogeny class, by a theorem of Kenku. □ 

This modular approach allows us, as a by-product, obtain a new bound for 
the cardinality of Y(Zs). 
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( 82 ) 


Theorem 5 . 0 . 3 . Let S contain all prime divisors of 26. Let 

Mi = JJ p 2+d r 
pes 

where d 2 = 6, ds = 3, and d p = 0 for p > 5. Let Xq(Mi) be the modular curve 
of level Jo (Mi), and let go(Mi) be its genus. Then, 

\Y(Z S )\ < 8 x |Aut Q (Y)| x <?o(Mi). (83) 

Note that go{M{) < M\. 

Proof. The dimension of space of cusp forms of weight two on Xo(Mi) is equal 
to the genus of Xq{M\). For each rational Hecke eigenform of weight two, there 
is at most eight isomorphism classes of elliptic curves by Kenku’s theorem. For 
each elliptic curve, there is at most |AutQ(V)| elements of F(Zg). Thus we 
obtain the claimed bound. □ 

Remark 5 . 0 . 1 . Above bound is clearly worse than previous ones, as go{M{) is 
roughly M\. Nevertheless, the proof of above bound uses modularity theorem 
as its key ingredient, and it has little to do with estimation of number of points 
on the curve y 2 = x 3 + n. 

In order to compute Y(Zg) in practice, one has to first tabulate elliptic 
curves. By a tabulation of elliptic curves, we shall mean the table of isomorphism 
classes of elliptic curves, represented in a Weierstrass equation, ordered by their 
conductors. In a sense this step of tabulation is a pre-computation, which only 
depends on the discriminant of h(x,y), and the table can be used again and 
again. 


6 Numerical examples 

We give numerical examples in this section. As explained before, the crucial step 
in working out a numerical example is to tabulate elliptic curves with specified 
conductor. We avoid this step by relying on Cremona’s Elliptic Curve Database. 
In particular, we assume throughout: 

Hypothesis. Cremona’s Elliptic Curve Datebase is conrplete(i,e., no curve 
is omitted) up to conductor 350,000. 

We will compute Y(Zs) based on Cremona’s database. The completeness of 
the list of the solutions depends on the truth of the hypothesis. More specifically, 
we need a complete list of elliptic curves which has good reduction outside of 
26 and S. 

We give some justification for introducing the above hypothesis. If some 
elliptic curves with conductor less than 350, 000 turn out to be omitted in the 
Cremona’s table, then it is possible that these new elliptic curves give rise to 
new solutions, which is not listed in the present paper. Such corrections can 
be always made upon each discovery, if any, of omitted elliptic curves. On the 
other hand, it is clearly beyond the scope of current work to check that the 
computer codes which were used to generate Cremona’s table contain no bugs. 


17 



6.1 Implementation and perfomance 

In this subsection, we explain how we implement the algorithm into a computer 
algebra package, and its performance. What we do not compute is the necessary 
table of elliptic curves. We assume that a list of elliptic curves of specified 
conductor is already available. 

In order to faithfully follow the proof of finiteness of Y(Zg), we should 
compute a set of elliptic curves which contains k(Y(Zs)), and compute n~ 1 {E) 
for each curve E in the set. For example, one might choose to compute the set 
of elliptic curves whose conductor divides the worst conductor of X t . However, 
this is computationally inefficients because of the following reason. When we 
replace E by its quadratic twist, Je(x,y 7 h, E) is replaced by its multiple. It 
follows that one has to solve the same polynomial 2 S+1 times when s is the 
cardinality of S. It turns out that working with X' t we avoid this problem of 
repeating J e (x,y, h, E). 

The following proposition tells us why it is possible to compute Y(Zg) using 
X' t instead of X t . 

Proposition 6.1.1. Let E be an elliptic curve. There exists t € Y(Zg) such 
that X t is isomorphic to a quadratic twist of E if and only if Jq(x, y, h , E) has 
a rational solution. 

Proof. Recall that we proved that Jq(x, y , h, E ), (l56l) . and (15711 have a common 
solution if and only if E is isomorphic to X t itself, without a quadratic twist. 
Hence, only if part is obvious, and we are left to prove if part of the proposition. 
It suffices to show that given a solution of Je(x,y,h, E), one can find a quadratic 
twist E' of E such that jQ(x,y, h,E'), (l56l) . and (1571) have a common solution. 

In order to see that finding such a quadratic twist E' is possible, recall that 
J 24 (x. y, h, E) was defined by 

J- 2 i{x, y, h, E) := A 12 (c 6 (E) 2 c 4 (x, y, h) 3 - c 4 (E) 3 c 6 (x, y, h ) 2 ) . (84) 

Let E is given by the equation 

E: y 2 = x 3 + a 4 x + a§ (85) 

then let E' be the quadratic twist 

E ': y 2 = x 3 + a 4 r 2 x + a^r 3 ( 86 ) 

of E by r. Then from the formula for J 2 4 {x, y,h, E), together with the relations 
m and m, tells us that 


J 24 (x, y, h, E') = r 6 ■ J 24 {x, y, h, E). (87) 

In particular, from factorisation (l59l) . it follows that 

J 6 {x, y, h, E') = r 6 • J 6 (x, y, h, E'). ( 88 ) 
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So what remains is to find a value of r, for a given zero (xq, yo ) of Jq{x, y , h, E') 
(or jQ(x,y,h,E)), for which 


C 4 {xo,yo,h ) = A 4 c 4 (U') (89) 

c 6 {x 0 ,yo, h ) = X 6 c 6 (E') (90) 


have a rational solution in A. Above equations are nothing but the equations 
(1551) and (1571) written for E'. Solving for A, we get 


A 2 _ c e (x 0 ,y 0 ,h) c 4 (E') 

c 4 (x 0 ,y 0 ,h) c 6 (E') 

_ c e {xo,y 0 ,h ) c 4 (E) 

C4(x 0 ,yo,h) cq(E) 


(91) 

(92) 


so one can take 


c 6 (x 0l y 0 ,h) c 4 (-E) 

c 4 {xo,yo, h) c 6 (E) 


(93) 


in order to render (l89l) and (l90l) to have solution A = ±1. The assertion of the 
proposition is proved. □ 

Based on Pronosition lG.l.ll we proceed as follows in order to compute Y(Zs). 
Suppose we are given h(x,y) and S. From this, one can compute the list of 
conductors of X[ for t £ Y( Zg), with negligible computational cost. Using 
the Cremona’s Elliptic Curve Database, we get a sequence of elliptic curves 
Ui, £ 2 , • ■ • ,Ek, for some finite fc, where each A, is given in a Weierstrass form. 
Now we compute rational solutions of Je(x,y,h,Ei) for each U,;. If there is a 
solution, say Xt and yt, then we proceed to compute h(xt,yt ) and double check 
that the result is correct. 

As we observed before, the coefficients of Je(x,y, h,E) has a large common 
factor when we start from an elliptic curve with a± = 0,2 = 0,3 = 0, we take 


J' 6 (x,y,h,E) = 2 22 -3 3 • J 6 (x,y,h,E) 


(94) 


whose coefficients are still integral. 

We give three examples of h, and present corresponding Jq(x , y, h, E), for a 
generic elliptic curve 


E: y 2 = x 3 + CL 4 X + ae- 


(95) 


Take 


hi(x,y) =x 2 y-xy 2 = x(x - y)y (96) 

which corresponds to the unit equation. In this case, 

J'&{x, y, hi,E) = (x- 2 y) 2 (x + y) 2 ( 2 x - y) 2 a 3 + 27{x 2 - xy + y 2 fal (97) 
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which has a particularly nice factorisation. In this case, the solutions (2,1), 
( 1 ,- 1 ), and ( 1 , 2 ), for h\{x,y) = ± 2 ° is clearly visible from the factors of 
coefficients of 04 , corresponding to curves with a 9 = 0. Asa less trivial example, 
we take 


: x 3 — x 2 — 21345a: — 1190943 


(98) 


which is the minimal Weierstrass equation for elliptic curve ”960e6” in Cre¬ 
mona’s database. We make change of variables x M- x + 1/3 


E: y 2 = x 6 - 


64036 


32347568 

27 


(99) 


for which 


J^x,y,h u E) (100) 

= - 64(3a; - 128y)(3a; + 125y)(125x - 128y)(125a; + 3y)(128a; - 125y)(128a; - 3 y) 

( 101 ) 

and indeed (3,128) belongs to Y(Zg) for S = {2,3,5}, because 128 = 2 7 , and 
128 — 3 = 125 = 5 3 . The other five factors corresponds to orbits of AutQ(Y). 

As an example for which k~ 1 (E) has no rational point, we take 

£ 960e5 : y 2 = x 3 - x 2 - 18465a; + 971937 (102) 

which is the curve ”960e5”, which is just next to the previous one. After routine 
change of variables x n- x + 1/3, we obtains 

Je{x,y,hi, E) (103) 

= - 64(25a; 2 - 13874au/ + 25y 2 )(25x 2 + 13824 xy - 13824a/ 2 )(13824a; 2 - 13824a :y - 25 y 2 ) 

(104) 

which is product of three irreducible quadratic polynomials. It follows that 
there is no solution t GY (Z 5 ) with S = {2, 3, 5} for which X[ is isomorphic to 
E 9 QQe5- 

Let us consider a different cubic form 


hi{x, y) = x 2 y + 7a / 3 = (a ; 2 + 7 y 2 )y 
which corresponds to Ramanujan-Nagell equation. In this case, 

4(a :,y,h 2 ,E) = 2 2 7 4 y 2 (9a; 2 + 7 y 2 ) 2 a\ - 3 3 7 3 (3a ; 2 - 7y 2 ) 3 a 2 e 
has again a nice factorisation. Take 

£21 oei : y 2 + xy = x 3 + 210a’ + 900 


which is isomorphic to 


y 2 = x 3 


10079 762481 

48 2 + 864 


(105) 

(106) 

(107) 


(108) 
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for which 


4(x,y,h 2 ,E) (109) 

= 1^24^ 45X ~ 47y )( 45a: + 47y)(2048a; 2 - 14805xy + 113561y 2 )(2048a; 2 + U805xy + 113561y 2 ). 

( 110 ) 

The solution (47,45) corresponds to 

/i 2 (47, 45) = 737280 (111) 

= 2 14 • 3 2 • 5 (112) 

which is an element of Y(Zs) for S = {2,3,5}. 

Finally, take 


h 3 (x,y) = x 3 - x 2 y - Axy 2 - y 3 (113) 

whose discriminant is 13 2 . In this case, 

Jq(x, y , h 3 , E) = 13 2 (5x 3 + 21a; 2 y + 6 xy 2 — 5y 3 ) 2 a 3 + 3 3 13 3 (a; 2 + xy + y 2 ) 3 Og 

(114) 

which shows that the coefficients of a 3 in general have irreducible factor of degree 
three. Of course, the squares and cubes in the coefficients are expected, as we 
have defined Je{x,y,h, E) by (158} and (l59l) . 

In practice, the coefficients 04 and a§ are very large, so it is difficult to factor 
Jg(a;,y, h,E) manually. Nonetheless, using suitable computer algebra packages, 
one can factor such polynomials rather quickly. The author’s experience shows 
that SageMathCloud is able to factor roughly 200 polynomials per second. 

Take h(x, y) = hi(x, y ) as above, and take S = {2,7,11,13}. Then we need 
a table of elliptic curves whose conductor divides 2® • 7 • 11 • 13 = 256256. Since 
256256 < 350000 we may use Cremona’s Elliptic Curve Database, from which 
we get 940 such curves. From them, we get 51 solutions, as we display in Tabled} 
It took 4.13 seconds in CPU time to generate Tabled] A typical box in the table 
looks like 


h(x t ,y t ) 

Factorisation of h(xt,yt) 
Cremona’s label 
Factorisation of the conductor 


(115) 


with five items listed vertically. 

Let us take S = {2,3, 5, 7} for h 2 (x, y). In this case, we need elliptic curves 
of conductor dividing 2® • 3 • 5 • 7 2 = 188160, which is smaller than 350000. The 
Cremona’s table contains 4568 such curves, from which we find 33 solutions. It 
took 20.09 seconds in CPU time to generate Tabled] 
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Let us take S = {2,5,13} for h 3 (x,y). In this case, we need elliptic curves 
of conductor dividing 2 8 • 5 • 13 2 = 216320, which is smaller than 350000. The 
Cremona’s table contains 976 such curves, from which we find 35 solutions. It 
took 4.20 seconds in CPU time to generate Table [T3l Since we wrote a computer 
code which compute the solutions with y t ^ 0, the trivial solution /i(l,0) = 1 
is omitted from the table . 

Remark 6.1.1. Repeated 6 ’s on the exponent of 2 in the conductors appearing 
in Table |T] can be perhaps predicted by considering connected components in 
the Neron model. On the other hand, as the 2 divides the discriminant of 
h- 2 (x, y), various exponents of 2 appear in Tabled] Although we ignored further 
analysis of conductors of X' t at the primes dividing 2 <5, such an analysis might 
help practical computations. 

6.2 Statistics for x(x — y)y 

In this section, we fix 


h{x,y) = x{x - y)y (116) 

and vary S' in a few directions. Let us begin with the case when 

S = {2,p} (117) 

consists of 2 and one more prime p. Cremona’s database allows us to compute 

Y{Z S ) if 

2® p< 350000 (118) 

or p < 1367. It follows that except for p = 5,7,17, 31,257, we have 

Y(Z S ) = Y( Z {2} ) = {(2 : 1 ), (1 : - 1 ), (1 : 2 )}. (119) 

In fact, this case is less interesting since computation of Y(Zs) reduces to 

2 m — p n = ±1 ( 120 ) 


which is a special case of Catalan’s equation. Since Catalan’s conjecture is 
known, solutions of above equation necessarily satisfies n = 1, and Y (Z,g) has 
more than three elements if and only if p is a prime of the form 2 m ± 1. Thus, 
we are merely verifying Catalan’s conjecture. 

As another example, take 


5 = {2,3,p} 


consists of 2, 3 and another prime p > 3. Since 


350000 
2® • 3 


< 456 


( 121 ) 


( 122 ) 
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one can use Cremona’s database for primes p up to 449, which is the 87th prime 
number. In the next table we display 

p : m 


as p varies among the primes from 5 to 449, and Y(Z $) has 3 + 6 m solutions. 


5 : 16 

7 : 12 

11 : 9 

13 : 8 

17 : 8 

19 : 7 

23 : 6 

29 : 5 

31 : 5 

37 : 5 

41 : 5 

43 : 5 

47 : 5 

53 : 4 

59 : 4 

61 : 5 

67 : 4 

71 : 4 

73 : 6 

79 : 4 

83 : 4 

89 : 4 

97 : 5 

101 : 4 

103 : 3 

107 : 4 

109 : 4 

113 : 4 

127 : 4 

131 : 4 

137 : 4 

139 : 4 

149 : 3 

151 : 3 

157 : 3 

163 : 4 

167 : 3 

173 : 3 

179 : 4 

181 : 3 

191 : 4 

193 : 4 

197 : 3 

199 : 3 

211 : 4 

223 : 3 

227 : 4 

229 : 4 

233 : 3 

239 : 4 

241 : 4 

251 : 4 

257 : 4 

263 : 3 

269 : 4 

271 : 3 

277 : 3 

281 : 3 

283 : 4 

293 : 3 

307 : 4 

311 : 3 

313 : 3 

317 : 3 

331 : 3 

337 : 4 

347 : 3 

349 : 3 

353 : 3 

359 : 3 

367 : 3 

373 : 3 

379 : 3 

383 : 4 

389 : 3 

397 : 3 

401 : 3 

409 : 3 

419 : 3 

421 : 3 

431 : 5 

433 : 4 

439 : 3 

443 : 3 

449 : 3 

// 

// 

// 

// 

// 


Note that for S = { 2 ,3}, Y(Z$) has 21 elements, so in particular m = 3. Above 
table indicates that m stabilises around 3, with a notable exception for p = 431. 
It is mainly due to a rather surprising identity 431 = 2 9 — 3 4 . 

Now we take 


S = {2,3,5,p} (123) 


where p is a prime number which does not exceed 89. Note that k/ 2 , 3 , 5 } contains 
99 = 3 + 6 x 16 elements, so a trivial lower bound for to in this case is 16. Using 
our algorithm, we obtain the following table. 


7 : 

62 

11 

46 

13 

44 

17 

37 

19 

37 

23 

35 

29 

31 

31 

30 

37 

30 

41 

30 

43 : 28 

47 

26 

53 

28 

59 

25 

61 

26 

67 

26 

71 

25 

73 

25 

79 

25 

83 

25 

89 

23 

// 


It indicates that Y(Zs) steadily decreases as p increases, although it is un¬ 
clear whether it will reach m = 16 at some point. Note an exceptional increment 
at p = 53, for which we record the solutions and associated Cremona label in 
Table [3] In fact, the number of possible elliptic curves tends to decrease as p 
increases. For example, if p = 7, there are 1688 curves, while the corresponding 
number is 1080 for p = 87. 

One wonders weather one can improve the bound on the cardinality of 
Y(Zs), as S varies among certain subsets of prime numbers with fixed car¬ 
dinality, such as S = So U {p} with varying p. 

7 Comparison with the work of Tzanakis and de 
Weger. 

There had been an attempt to explicitly solve Thue-Mahler equation by Tzanakis 
and de Weger, based on linear forms in logarithms. Theoretical foundations for 
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two approaches are quite different, and in this section we compare the two from 
practical point of view. 

In our approach, the computation of k _ 1 (£ i ) for a given E is easy. There is a 
formula for Jq{x , y, h, E) to which we plug in the coefficients of E , and the ratio¬ 
nal solutions of Jq(x , y 1 h , E) can be found quickly. On the other hand, finding 
all possible candidates of E is difficult, although the modularity of elliptic curves 
significantly facilitates it. The upshot is that in our approach, the modularity of 
elliptic curves reduces the determination of all possible elliptic curves to a single 
linear algebra problem on the space of modular forms. One can model the space 
of modular forms using the space of modular symbols. As Cremona explains 
in Section 2 of his article j§,, computation of modular symbols of a given level 
can be done rather quickly. On this Q-vector space of modular symbols, whose 
dimension is quite large, one has to find all one-dimensional Hecke invariant 
rational subspaces, and compute sufficiently many Hecke eigenvalues in order to 
find approximated C 4 and cq invariants of a curve in the corresponding isogeny 
class. Given an isogeny class, it is not so difficult to determine all isomorphism 
classes which belong to it. To summarise, this linear algebra problem on a huge 
space of modular symbols seems to lie at the bottleneck of our process. 

The approach of Tzanakis and de Weger is, as explained in the introduction 
of US]: consists of three steps. The first is to obtain large bounds from estimation 
of linear forms in logarithms of possibly irrational algebraic numbers. The 
second is to reduce the large bounds using diophantine approximation. The last 
step is to search for solutions below the bound, not by brutal force, but by using 
an algorithm to search for lattice points on a given sphere, a sieving process, 
and enumeration of possibilities. The authors remarks that the third process 
might well be the computational bottleneck for their process. They worked out 
two following concrete examples 

a: 3 - 23 x 2 y + 5 xy 2 + 24 y 3 = ±2 ei 3 e2 5 e3 7 64 , 5 = 5 2 • 44621 (124) 

x 3 - 3 xy 2 - y 3 = ±3 ei 17 62 19 63 , S = 3 4 (125) 

using their method. In order to solve above equation using our method, we 
need elliptic curves with conductor dividing 2 8 • 3 ■ 5 2 • 7 • 44621 2 > 5 x 10 13 , and 
2 8 ■ 3 5 • 17 • 19 > 2 x 10', which are not provided by Cremona’s database. 

We observe that the tools of our method has little to do with those of 
Tzanakis and de Weger. The computational bottlenecks of two approaches 
are different: ours is in linear algebra while theirs (seems to) be in geometry of 
numbers. As we vary h(x,y), we observe another difference. The computations 
we need to carry are mostly insensitive to h(x, y) except for the discriminant and 
S. Once the database of elliptic curves is established one can use the same data 
for a different h(x, y). The computations of Tzanakis and de Weger depends on 
the specific S'-unit equation which is sensitive to a chosen zero of h{x,y). 
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8 Generalised Ramanujan-Nagell equation 

The goal of current section is twofold. Firstly, we shall consider a special form of 
h(x , y) and determine Y (Zg), from which we deduce the complete set of solutions 
of certain generalised Ramanujan-Nagell equations. Secondly, we shall analyse 
the statistical behaviour as we vary S as h{x , y) remains fixed. 

The equation 


x 2 + 7 = T (126) 

for integers x and n , is often called the Ramanujan-Nagell equation in the 
literature. One can relate it to Thue-Mahler equation, since if we take 

h(x,y) = (x 2 +dy 2 )y (127) 

then the solution (x, n) = (xo,no) of the Ramanujan-Nagell equation leads to 
the point (x,y) = (xq, 1) £ Y(Zg) for d = 7 and S = {2}. Conversely we can 
recover the solutions of the Ramanujan-Nagell equation by computing Y(Zg). 
Thus one may consider Thue-Mahler equations for (1 1271) as a generalisation of 
Ramanujan-Nagell equation. 

In Table [H we display the elements of Y(Zg) for d = 7 and S = {2,3, 5, 7}, 
from which we conclude 


x 2 + 7 = 2 ei 3 e2 5 e3 7 e4 (128) 

with positive x has seven solutions corresponding to x = 181, 21,11, 7, 5, 3 and 
1. In Table 0 we find that 


x 2 + 7 = 2 ei 7 e2 ll e3 13 e4 (129) 

in positive x has fourteen solutions corresponding to x = 273, 181, 75, 53, 35, 
31, 21, 13, 11, 9, 7, 5, 3, and 1. 

For d = —7, and S = {2, 5, 7,11}, Table fl2l shows that we have particularly 
few solutions. Indeed Y( Zg) has five elements, among which three elements 
satisfy y = 1. In particular, 


x 2 — 7 = 2 ei 5 e2 7 e3 ll 64 (130) 

has only one solution x = 3 among positive integers. 

As we vary S for a fixed d, such that —d is not a square, we observe a pattern 
which we describe now. The pattern seems to persist for any such d, but let us 
take d = 1 for clarity. In particular, we consider 

h(x,y) = {x 2 + y 2 )y (131) 

in the rest of the present section. 

Take S = {2,p}, for a prime p > 3. Let the cardinality of Y(Z{ 2 ,p}) to be m. 
We shall divide it into two cases, depending on whether or not — 1 is a quadratic 
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residue modulo p, and compare the variation of in. Note that Y (Z{ 2 }) has three 
elements corresponding to x = 0,1 and —1. The table below lists 


p : m 


in the range of p for which —1 is a quadratic residue, and p < 113. 


5 : 15 

13 : 9 

17 : 9 

29 : 7 

37 : 5 

41 : 7 

53 : 5 

61 : 5 

73 : 5 

89 : 5 

97 : 5 

101 : 5 

109 : 3 

113 : 7 


(132) 


On the other hand, we observe that 

Y(Z {2>p} )=Y(Z {2} ) (133) 

if —1 is quadratic non-residue modulo p and p < 113. 

The different behaviour of the cardinality of Y(Zg) continues when we take 
S = {2,3,p}. We can numerically verify that 

y(Z {2)3 , p} ) = Y(Z {2} ) (134) 

if —1 is quadratic non-residue modulo p , and p < 455. In contrast, if —1 is a 
quadratic residue modulo p and p < 455, then Y( Z{ 2 , 3 , P }) is strictly larger than 
Y( Z{ 2 }) all the time. For instance, for p = 449, 

(13 2 + 27 2 ) • 27 = 24246 = 3 • 3 3 • 449 (135) 

is associated to the curve ”172416ol”, and similarly for p = 433, 

(17 2 + 12 2 ) • 12 = 5196 = 2 4 • 3 • 433 (136) 

is associated to the curve ”20784e2”. The numerical data suggests that the 
cardinality of Y (Z^ 23p y) is exactly five for most of p. 

Now we take S = {2, 5 ,p}. Note that Y(Z{ 2 , 5 }) consists of fifteen elements. 
Based on the previous observation, one might conjecture that Y(Zr 2 5p j) con¬ 
sists of fifteen elements if —1 is quadratic non-residue module p. However, there 
is a counterexample for p = 139. Indeed, one finds 

(29 2 + 278 2 ) • 278 = 2 • 5 7 • 139 (137) 

which is associated to the curve ”11120e2”. Among the primes 5 < p < 271 for 
which —1 is a quadratic non-residue, we verify 

Y(Z {2 , 5tP} )=Y(Z {2M ) (138) 

holds exceptp = 7,11,19,31, 79,139,191. In these exceptional cases, Y (Z{ 2i5 , p }) — 
y(Zr 2; 5 }) contains four elements when p = 7,11, and two elements in the re¬ 
maining five cases. In contrast, if we take S = {2, 5 ,p} for a prime p for which 
— 1 is a quadratic residue, then Y(Z^ 2 5 p y) is strictly larger than y(Zr 2 5 \), in 
the range 5 < p < 271. The smallest cardinality of Y( Z{ 2i5iP j) is seventeen, 
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which happens precisely for p = 241. Up to sign of Xt there is a unique element 
of Y(Z{ 2,5,241}) which does not belong to Y(Z{ 2 , 5 }), which is 

(15 2 + 4 2 )4 = 964 = 2 2 • 241 (139) 

associated to the curve ”15424d2”. 

We conclude by formulating a precise question based on our observation. 
Let So be a fixed set of primes. For a positive real number X and i £ {1,3}, let 
ni(X) be the number of primes p less than X which are congruent to i modulo 
4. Define 


Ai(S'o) = liminf 


X—>oo TTi(X) \ ^ 

' p<X,p=l(mod 4) 




and 


(140) 


^3 (So) = lirnsup ——— 
A'->oo 7^3 (A J 


X! \ Y (ZSoU{ P })\ 

^p<-X',p=3(mod 4) 


(141) 


One might speculate that ^4i({2}) = 5 and Aa({2})) = 3. To be on a conserva¬ 
tive side, one might ask whether 


A^So) > A 3 (S 0 ) 

holds. The author is not able to show 

bll(So) > ^3(So) 


(142) 


(143) 


holds for any particular So- 
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Table 1: (a, b, c, d) = (0,1, -1,0), S = 1, S = (2, 7,11,13} 


(-343,169) 
29679104 

2 9 ■ 7 3 ■ 13 2 
5824bd2 

2 6 ■ 7 ■ 13 

(-169,7) 

208208 

2 4 ■ 7 - 11■13 2 
64064a2 

2® ■ 7- 11 ■ 13 

(-169,343) 

29679104 

2 9 ■ 7 3 ■ 13 2 
5824bd2 

2® ■ 7 ■ 13 

(-121, 7) 

108416 

2 7 - 7 - ll 2 

4928h2 

2® ■ 7 - 11 

(-64, 13) 

64064 

2® ■ 7- 11 ■ 13 
64064k2 

2® ■ 7 - 11 ■ 13 

(-13,1) 

(-13,64) 

(-H.2) 

(-7,1) 

(-7,4) 

182 

64064 

286 

56 

308 

2-7-13 

2® ■ 7- 11 ■ 13 

2-11-13 

2 3 - 7 

2 2 ■ 7- 11 

5824b2 

64064k2 

9152bd2 

448b2 

4928m2 

2 6 ■ 7 ■ 13 

2® ■ 7 - 11 ■ 13 

2® - 11 • 13 

2® - 7 

2® ■ 7- 11 

(-7,121) 

(-7,169) 

(-4,7) 

(-2,11) 

(-1,1) 

108416 

208208 

308 

286 

2 

2 7 - 7 ■ ll 2 

2 4 ■ 7 - 11■13 2 

2 2 ■ 7- 11 

2-11-13 

2 

4928h2 

64064a2 

4928m2 

9152bd2 

256c2 

2 6 ■ 7- 11 

2® ■ 7 - 11 ■ 13 

2® ■ 7- 11 

2® ■ 11 ■ 13 

2® 

(-1,7) 

(-1,13) 

(1,2) 

(1,8) 

(1,14) 

56 

182 

-2 

-56 

-182 

2 3 • 7 

2-7-13 

-1 ■ 2 

-1 ■ 2 3 - 7 

-1-2-7-13 

448b2 

5824b2 

256c2 

448b2 

5824b2 

2 6 • 7 

2® - 7 - 13 

2 s 

2® ■ 7 

2® ■ 7 ■ 13 

(2,1) 

(2, 13) 

(4, 11) 

(7, 8) 

(7,11) 

2 

-286 

-308 

-56 

-308 

2 

-1-2-11-13 

-1 • 2 2 ■ 7 - 11 

-1 ■ 2 3 - 7 

-1 • 2 2 ■ 7 - 11 

256c2 

9152bd2 

4928m2 

448b2 

4928m2 

2® 

2® ■ 11 - 13 

2® ■ 7- 11 

2® ■ 7 

2® ■ 7 - 11 

(7,128) 

(7,176) 

(8,1) 

(8,7) 

(H,4) 

-108416 

-208208 

56 

56 

308 

-1 • 2 7 • 7- ll 2 

-1 - 2 4 ■ 7 ■ 11 ■ 13 2 

2 3 - 7 

2 3 ■ 7 

2 2 ■ 7 - 11 

4928h2 

64064a2 

448b2 

448b2 

4928m2 

2 6 ■ 7 - 11 

2® ■ 7 - 11 ■ 13 

2® - 7 

2® - 7 

2® ■ 7- 11 

(11,7) 

(11.13) 

(13,2) 

(13,11) 

(13,14) 

308 

-286 

286 

286 

-182 

2 2 ■ 7- 11 

-1-2-11-13 

2-11-13 

2-11-13 

-1-2-7-13 

4928m2 

9152bd2 

9152bd2 

9152bd2 

5824b2 

2 6 ■ 7 - 11 

2® ■ 11 - 13 

2® ■ 11 • 13 

2® ■ 11 ■ 13 

2® ■ 7 - 13 

(13, 77) 

(14, 1) 

(14,13) 

(64, 77) 

(77,13) 

-64064 

182 

182 

-64064 

64064 

-1 ■ 2 6 ■ 7- 11 ■ 13 

2-7-13 

2-7-13 

-1 -2® • 7 - 11 • 13 

2® ■ 7■11 - 13 

64064k2 

5824b2 

5824b2 

64064k2 

64064k2 

2® ■ 7 ■ 11 ■ 13 

2® ■ 7 - 13 

2® ■ 7 - 13 

2® - 7- 11•13 

2® ■7 ■ 11 - 13 

(77, 64) 

(121,128) 

(128, 7) 

(128,121) 

(169,176) 

64064 

-108416 

108416 

108416 

-208208 

2® ■ 7 ■ 11 ■ 13 

-1 ■ 2 7 - 7 - ll 2 

2 7 -7 - ll 2 

2 7 -7 - ll 2 

-1 ■ 2 4 ■ 7- 11 - 13 2 


4928h2 

4928h2 

4928h2 

64064a2 

2® ■ 7 ■ 11 ■ 13 

2® - 7 - 11 

2® ■ 7 - 11 

2® ■ 7 ■ 11 

2® - 7■11 - 13 

(169,512) 

(176, 7) 

(176,169) 

(343, 512) 

(512,169) 

-29679104 

208208 

208208 

-29679104 

29679104 

-1 ■ 2 9 • 7 3 • 13 2 

2 4 - 7 - 11 - 13 2 

2 4 • 7 ■ 11 ■ 13 2 

-1 - 2 9 ■ 7 3 - 13 2 

2 9 ■ 7 3 - 13 2 

5824bd2 

64064a2 

64064a2 

5824bd2 

5824bd2 

2® ■ 7 ■ 13 

2® ■ 7 - 11 ■ 13 

2®■7 ■ 11 - 13 

2® ■ 7 ■ 13 

2® ■ 7 - 13 

(512, 343) 
29679104 

2 9 ■ 7 3 ■ 13 2 

empty 

empty 

empty 

empty 

5824bd2 

2® ■ 7 ■ 13 
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Table 2: (a, 6, c, d) = (0,1, -1,0), 6 = 1, S = {2,3,431} 


(-431,1) 

(-431,81) 

(-81,431) 

(-8,1) 

(-3,1) 

(-2,1) 

186192 

17874432 

17874432 

72 

12 

6 

2 4 ■ 3 3 ■ 431 

2 9 ■ 3 4 - 431 

2 9 ■ 3 4 ■ 431 

2 3 ■ 3 2 

2 2 - 3 

2 ■ 3 

82752bd2 

82752t2 

82752t2 

192c3 

192d2 

192b2 

2® • 3•431 

2®- 3 ■ 431 

2® • 3•431 

2® • 3 

2® - 3 

2® • 3 

(-1,1) 

(-1,2) 

(-1,3) 

(-1,8) 

(-1,431) 

(1,2) 

2 

6 

12 

72 

186192 

-2 

2 

2 ■ 3 

2 2 • 3 

2 3 ■ 3 2 

2 4 ■ 3 3 • 431 

-1 ■ 2 

256c2 

192b2 

192d2 

192c3 

82752bd2 

256c2 

2® 

2® • 3 

2® • 3 

2® • 3 

2® ■ 3 ■ 431 

2 s 

(1,3) 

(1,4) 

(1,9) 

(1,432) 

(2,1) 

(2,3) 

-6 

-12 

-72 

-186192 

2 

-6 

-1-2-3 

-1 ■ 2 2 ■ 3 

-1 ■ 2 3 - 3 2 

-1 ■ 2 4 ■ 3 3 - 431 

2 

-1-2-3 

192b2 

192d2 

192c3 

82752bd2 

256c2 

192b2 

2® • 3 

2® • 3 

2® • 3 

2®- 3 - 431 

2 s 

2® - 3 

(3, 1) 

(3,2) 

(3,4) 

(4,1) 

(4,3) 

(8,9) 

6 

6 

-12 

12 

12 

-72 

2 • 3 

2 • 3 

-1 ■ 2 2 ■ 3 

2 2 ■ 3 

2 2 - 3 

-1 • 2 3 ■ 3 2 

192b2 

192b2 

192d2 

192d2 

192d2 

192c3 

2® ■ 3 

2® • 3 

2® • 3 

2® ■ 3 

2® ■ 3 

2® - 3 

(9, 1) 

(9,8) 

(81,512) 

(431, 432) 

(431,512) 

(432,1) 

72 

72 

-17874432 

-186192 

-17874432 

186192 

2 3 • 3 2 

2 3 • 3 2 

-1 ■ 2 9 ■ 3 4 -431 

-1 ■ 2 4 ■ 3 3 - 431 

-1 - 2 9 - 3 4 ■ 431 

2 4 ■ 3 3 - 431 

192c3 

192c3 

82752t2 

82752bd2 

82752t2 

82752bd2 

2® ■ 3 

2® • 3 

2®- 3 - 431 

2®- 3 - 431 

2® ■ 3■431 

2®■3 ■ 431 

(432, 431) 

(512,81) 

(512,431) 




186192 

17874432 

17874432 




2 4 ■ 3 3 ■ 431 

2 9 ■ 3 4 - 431 

2 9 ■ 3 4 - 431 

empty 

empty 

empty 

82752bd2 

82752t2 

82752t2 




2®■3 ■ 431 

2®- 3 ■ 431 

2® • 3■431 
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Table 3: (a, 6, c, d) = (0,1, -1,0), 6 = 1, S = (2,3,5,53} 


(-6561,64) 

50880cd2 

(-3072, 53) 
50880bq2 

(-256,9) 

50880y2 

(-159, 1) 
50880bx2 

(-125,3) 

960e6 

(-81,25) 

50880bw2 

(-80,1) 

960g5 

(-75,53) 

50880dz2 

(-72,53) 

50880bv2 

(-64,6561) 

50880cd2 

(-53,1) 
10176f2 

(-53, 27) 
50880s2 

(-53,72) 

50880bv2 

(-53, 75) 
50880dz2 

(-53, 3072) 
50880bq2 

(-50,3) 

50880cw2 

(-48,5) 

50880bt2 

(-45,8) 

50880ea2 

(-27,5) 

960o2 

(-27, 53) 
50880s2 

(-25,2) 

960d2 

(-25,81) 

50880bw2 

(-24,1) 
960b3 

(-16,9) 

960g3 

(-15,1) 

960i2 

msxh 

(-9,16) 
960g3 

(-9,256) 

50880y2 

(-8,1) 

192c3 

(-8,45) 

50880ea2 

(-5,1) 

960f2 

(-5,3) 

960c2 

(-5,4) 

96012 

(-5,27) 

960o2 

(-5,48) 

50880bt2 

(-4,1) 

320b2 

(-4,5) 

96012 

(-3,1) 

192d2 

(-3,2) 

960n2 

(-3,5) 

960c2 

(-3, 50) 
50880cw2 

(-3,125) 
960e6 

(-2,1) 

192b2 

(-2,3) 

960n2 

(-2,25) 

960d2 

(-1,1) 

256c2 

(-1,2) 

192b2 

■H 

(-1,4) 

320b2 

.jggsM. 

WEMM 

(-1,9) 

960j2 

(-1,15) 

960i2 

(-1,24) 

960b3 

(-1.53) 

10176f2 

(-1,80) 

960g5 

(-1, 159) 
50880bx2 

(1,2) 

256c2 

(1,3) 

192b2 

(1,4) 

192d2 

(1,5) 

320b2 

(1,6) 

960f2 

(1,9) 

192c3 

(1,10) 

960j2 

(1,16) 

960i2 

(1,25) 

960b3 

(1,54) 

10176f2 

(1,81) 

960g5 

(1,160) 
50880bx2 

(2,1) 

256c2 

(2,3) 

192b2 

(2,5) 

960n2 

(2,27) 

960d2 

(3,1) 

192b2 

(3,2) 

192b2 

(3,4) 

192d2 

(3,5) 

960n2 

(3,8) 

960c2 

(3, 53) 
50880cw2 


(4,1) 

192d2 

(4,3) 

192d2 

(4,5) 

320b2 

(4,9) 

96012 

(5,1) 

320b2 

(5,2) 

960n2 

(5,3) 

960n2 

(5,4) 

320b2 

(5,6) 

960f2 


(5,9) 

96012 

(5,32) 

960o2 

(5, 53) 
50880M2 

(6,1) 

960f2 

(6,5) 

960f2 

[Ml* * W 

(8,5) 

960c2 

(8,9) 

192c3 

(8, 53) 
50880ea2 

(9,1) 

192c3 

(9,4) 

96012 

(9,5) 

96012 

(9,8) 

192c3 

(9,10) 

960j2 

(9, 25) 

960g3 

(9, 265) 
50880y2 

(10,1) 
960j2 

(10,9) 

960j2 

(15,16) 

960i2 

(16,1) 

960i2 

(16,15) 

960i2 


(24,25) 
960b3 

(25,1) 

960b3 

(25, 9) 

960g3 


(25,24) 
960b3 

(25, 27) 
960d2 

(25,106) 
50880bw2 

(27, 2) 

960d2 

(27, 25) 
960d2 


(27, 80) 
50880s2 

(32, 5) 

960o2 


(45, 53) 
50880ea2 

(48, 53) 
50880bt2 

(50, 53) 
50880cw2 

(53, 3) 
50880cw2 

(53, 5) 
50880bt2 

(53,8) 

50880ea2 

(53, 45) 
50880ea2 

(53,48) 

50880bt2 

(53, 50) 
50880cw2 

(53, 54) 
10176f2 

(53, 80) 
50880s2 

(53,125) 
50880bv2 

(53,128) 
50880dz2 

(53, 3125) 
50880bq2 

(54, 1) 
10176f2 

(54, 53) 
10176f2 

(64, 6625) 
50880cd2 

(72,125) 
50880bv2 

(75,128) 
50880dz2 


(80, 27) 
50880s2 


(106, 25) 
50880bw2 


(128, 53) 
50880dz2 


(256, 265) 
50880y2 


(6561,6625) 
50880cd2 


(80, 53) 
50880s2 


(106,81) 

50880bw2 


(128, 75) 
50880dz2 


(265, 9) 
50880y2 


(6625, 64) 
50880cd2 


(125,53) 

50880bv2 


(128. 125) 
960e6 


6 ) 


(6625, 6561) 
50880cd2 


(81, 1) 
960g5 


(125, 72) 
50880bv2 


(159, 160) 
50880bx2 


(3072, 3125) 
50880bq2 


(125,128) 
960e6 


(160, 1) 
50880bx2 


(3125, 53) 
50880bq2 


(81,106) 
50880bw2 


(128,3) 

960e6 


(160, 159) 
50880bx2 


(3125, 3072) 
50880bq2 










































































































































































Table 4: (a, b , c, d) = (0,1,0,7), 6 = -1 • 2 2 • 7, S' = {2,3,5,7} 


(-181,1) 

32768 

2 15 

3136y5 

2 6 • 7 2 

(-119,5) 

71680 

2 11 -5-7 
15680vl 

2 6 • 5 • 7 2 

(-47,45) 

737280 

2 14 • 3 2 •5 
210el 

2 • 3 • 5 • 7 

(-35,9) 

16128 

2 8 • 3 2 • 7 
294gl 

2 • 3 • 7 2 

(-31,3) 

3072 

2 10 -3 
9408pl 

2 6 • 3 • 7 2 

(-21,1) 

(-13,7) 

(-H.1) 

(-9,5) 

(-7,1) 

448 

3584 

128 

1280 

56 

2 6 • 7 

2 9 • 7 

2 7 

2 8 • 5 

2 3 • 7 

3136r3 

14al 

3136kl 

70al 

3136hl 

2 6 • 7 2 

2-7 

2 6 • 7 2 

2-5-7 

2 6 • 7 2 

(-7,3) 

(-7,5) 

(-5,1) 

(-3,1) 

(-1,1) 

336 

1120 

32 

16 

8 

2 4 • 3 • 7 

2 5 • 5 • 7 

2 5 

2 4 

2 3 

2352gl 

3920ol 

3136bbl 

112bl 

224al 

2 4 • 3 • 7 2 

2 4 •5 • 7 2 

2 6 • 7 2 

2 4 • 7 

2 5 • 7 

(-1.3) 

(0,1) 

(1,1) 

(1,3) 

(3,1) 

192 

7 

8 

192 

16 

2 6 -3 

7 

2 3 

2 6 -3 

2 4 

21a4 

12544a2 

224al 

21a4 

112bl 

3-7 

2 8 • 7 2 

2 5 • 7 

3 • 7 

2 4 • 7 

(5,1) 

(7,1) 

(7,3) 

(7,5) 

(9,5) 

32 

56 

336 

1120 

1280 

2 5 

2 3 • 7 

2 4 • 3 • 7 

2 5 • 5 • 7 

2 8 -5 

3136bbl 

3136hl 

2352gl 

3920ol 

70al 

2 6 • 7 2 

2 6 • 7 2 

2 4 • 3 • 7 2 

2 4 • 5 • 7 2 

2-5-7 

(11,1) 

(13,7) 

(21,1) 

(31,3) 

(35,9) 

128 

3584 

448 

3072 

16128 

2 7 

2 9 • 7 

2 6 • 7 

2 10 -3 

2 8 • 3 2 • 7 

3136kl 

14al 

3136r3 

9408pl 

294gl 

2 6 • 7 2 

2-7 

2 6 • 7 2 

2 6 • 3 • 7 2 

2 • 3 • 7 2 

(47,45) 

737280 

2 14 • 3 2 • 5 

(119,5) 

71680 

2 11 • 5 • 7 

(181,1) 

32768 

2 15 

empty 

empty 

210el 

2 • 3 • 5 • 7 

15680vl 

2 6 • 5 • 7 2 

3136y5 

2 6 • 7 2 
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Table 5: (a,b,c,d) = (0,1,0,7), 5 = -1 • 2 2 • 7, S = {2,7,11,13} 


(-273,1) 
74536 

2 3 • 7 • ll 3 
34496dcl 

2 6 • 7 2 ■ 11 

(-181,1) 

32768 

2 15 

3136y5 

2® ■ 7 2 

(-161,13) 
352352 

2 5 ■ 7■ll 2 ■ 13 
112112ol 

2 4 ■ 7 2 • 11 ■ 13 

(-87,91) 

5963776 

2 1 ® ■ 7 • 13 
182al 

2 ■ 7 - 13 

(-81,13) 
100672 

2® ■ ll 2 ■ 13 
112112bfl 

2 4 ■ 7 2 - 11 ■ 13 

(-75,1) 

5632 

2 9 - 11 
34496bfl 

2® - 7 2 ■ 11 

(-67, 49) 

(-57,11) 

(-53,1) 

(-49,13) 

(-35,1) 

(-31,1) 

1043504 

45056 

2816 

46592 

1232 

968 

2 4 ■ 7 2 ■ ll 3 

2 12 ■ 11 

2 8 ■ 11 

2 9 ■ 7■13 

2 4 ■ 7 - 11 

2 3 ■ ll 2 

1232dl 

34496cfl 

34496dol 

1274jl 

34496zl 

34496dj 1 

2 4 ■ 7■11 

2® ■ 7 2 ■ 11 

2® ■ 7 2 ■ 11 

2 ■ 7 2 ■ 13 

2® ■ 7 2 - 11 

2® - 7 2 ■ 11 

(-25, 7) 

(-21,1) 

(-19,7) 

(-15,13) 

(-13,1) 

(-13,7) 

6776 

448 

4928 

18304 

176 

3584 

2 3 ■ 7 ■ ll 2 

2® ■ 7 

2® ■ 7- 11 

2 7 ■ 11 ■ 13 

2 4 • 11 

2 9 - 7 

2464dl 

3136r3 

77cl 

2002cl 

34496csl 

14al 

2 5 • 7■11 

2® ■ 7 2 

7-11 

2 ■ 7 ■ 11 ■ 13 

2® - 7 2 - 11 

2 ■ 7 

(-H.1) 

(-9,1) 

(-7,1) 

(-7,2) 

(-7,11) 

(-7,13) 

128 

88 

56 

154 

9856 

16016 

2 7 

2 3 ■ 11 

2 3 ■ 7 

2 • 7 - 11 

2 7 ■7 - 11 

2 4 ■ 7 - 11 • 13 

3136kl 

34496bbl 

3136hl 

8624h2 

107811 

112112il 

2 6 ■ 7 2 

2® ■ 7 2 ■ 11 

2® ■ 7 2 

2 4 • 7 2 • 11 

2 - 7 2 ■ 11 

2 4 • 7 2 - 11 • 13 

(-5,1) 

(-3,1) 

(-3,4) 

(-3,7) 

(-2,1) 

(-1-1) 

32 

16 

484 

2464 

11 

8 

2 s 

2 4 

2 2 ■ ll 2 

2 5 ■ 7■11 

11 

2 3 

3136bbl 

112bl 

616e4 

1232el 

2464k1 

224al 

2 6 ■ 7 2 

2 4 • 7 

2 3 ■ 7 - 11 

2 4 ■ 7■11 

2 5 ■ 7- 11 

2 5 ■ 7 

(0,1) 

(1,1) 

(2,1) 

(3,1) 

(3,4) 

(3,7) 

7 

8 

11 

16 

484 

2464 

7 

2 3 

11 

2 4 

2 2 ■ ll 2 

2 5 ■ 7 - 11 

12544a2 

224al 

2464kl 

112bl 

616e4 

1232el 

2 8 ■ 7 2 

2 5 ■ 7 

2 5 ■ 7-11 

2 4 • 7 

2 3 ■ 7■11 

2 4 ■ 7■11 

(5,1) 

(7,1) 

(7,2) 

(7,11) 

(7,13) 

(9,1) 

32 

56 

154 

9856 

16016 

88 

2 5 

2 3 ■ 7 

2 ■ 7 - 11 

2 7 • 7■11 

2 4 - 7- 11 - 13 

2 3 - 11 

3136bbl 

3136hl 

8624h2 


11211211 

34496bbl 

2 6 ■ 7 2 

2® ■ 7 2 

2 4 • 7 2 ■ 11 

2 • 7 2 ■ 11 

2 4 ■ 7 2 - 11 ■ 13 

2® ■ 7 2 ■ 11 

(11,1) 

(13,1) 

(13,7) 

(15,13) 

(19,7) 

(21,1) 

128 

176 

3584 

18304 

4928 

448 

2 7 

2 4 • 11 

2 9 ■ 7 

2 7 ■ 11 ■ 13 

2® ■ 7 - 11 

2® ■ 7 

3136kl 

34496csl 

14al 

2002cl 

77cl 

3136r3 

2 6 ■ 7 2 

2® ■ 7 2 ■ 11 

2 • 7 

2 ■ 7 ■ 11 ■ 13 

7- 11 

2® ■ 7 2 

(25,7) 

(31,1) 

(35,1) 

(49,13) 

(53,1) 


6776 

968 

1232 

46592 

2816 


2 3 ■ 7 • ll 2 

2 3 • ll 2 

2 4 ■ 7 ■ 11 

2 9 ■ 7■13 

2 8 ■ 11 

2 12 ■ 11 

2464dl 

34496dj 1 

34496zl 

1274jl 

34496dol 

34496cfl 

2 s ■ 7 ■ 11 

2® ■ 7 2 ■ 11 

2® ■ 7 2 ■ 11 

2 • 7 2 ■ 13 

2® ■ 7 2 - 11 

2® ■ 7 2 ■ 11 

(67, 49) 

(75,1) 

(81,13) 

(87, 91) 

(161,13) 

(181,1) 

1043504 

5632 

100672 

5963776 

352352 

32768 

2 4 ■ 7 2 ■ ll 3 

2 9 • 11 

2® ■ ll 2 • 13 

2 1 ® • 7 • 13 

2 5 • 7■ll 2 ■ 13 

2 15 

1232dl 

34496bfl 

112112bfl 

182al 

112112ol 

3136y5 

2 4 ■ 7■11 

2® ■ 7 2 ■ 11 

2 4 ■ 7 2 • 11 ■ 13 

2-7-13 

2 4 ■ 7 2 - 11 ■ 13 

2® ■ 7 2 

(273,1) 
74536 

2 3 ■ 7 ■ ll 3 

empty 

empty 

empty 

empty 

empty 

34496dcl 

2® ■ 7 2 • 11 
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Table 6: (a, 6, c, d) = (0,1,0,2), S = -1 • 2 3 , S = {2,3,11} 


(-695, 8) 

(-241, 22) 

(-155,4) 

(-140,1) 

(-59,16) 

(-50,9) 

3865224 

1299078 

96228 

19602 

63888 

23958 

2 3 • 3 • 11 s 

2 - 3 10 ■ 11 

2 2 ■ 3 7 • 11 

2 • 3 4 ■ ll 2 

2 4 ■ 3 - ll 3 

2 • 3 2 - ll 3 

2112d4 

2112n4 

2112j2 

8448e2 

2112e4 

8448q2 

2 6 ■ 3 • 11 

2 6 - 3 - 11 

2® - 3 ■ 11 

2 8 • 3•11 

2® • 3 - 11 

2 8 • 3 - 11 

(-31,8) 

(-25,32) 

(-22, 1) 

(-19, 1) 

(-17,2) 

(-14,1) 

8712 

85536 

486 

363 

594 

198 

2 3 ■ 3 2 ■ ll 2 

2 5 ■ 3 s • 11 

2 - 3 5 

3 ■ ll 2 

2 ■ 3 3 ■ 11 

2 - 3 2 ■ 11 

2112r4 

66c2 

768d4 

4224s2 

2112b2 

8448f2 

2 6 ■ 3 • 11 

2-3-11 

2 8 • 3 

2 7 • 3•11 

2® • 3■11 

2 8 • 3■11 

(-13,8) 

(-8,1) 

(-7,4) 

(-7,6) 

(-5,1) 

(-5,2) 

2376 

66 

324 

726 

27 

66 

2 3 ■ 3 3 • 11 

2-3-11 

2 2 ■ 3 4 

2 ■ 3 ■ ll 2 

3 3 

2-3-11 

66a2 

8448c2 

24a6 

1056c4 

384g2 

211212 

2 • 3 • 11 

2 8 - 3 - 11 

2 3 ■ 3 

2® • 3•11 

2 7 - 3 

2® • 3 • 11 

(-4, 1) 

(-3,1) 

(-2,1) 

(-2,3) 

(-1,1) 

(-1,2) 

18 

11 

6 

66 

3 

18 

2 ■ 3 2 

11 

2 ■ 3 

2-3-11 

3 

2 • 3 2 

768c2 

1408d2 

768hl 

8448tl 

384b 1 

96b4 

2® • 3 

2 7 • 11 

2 8 ■ 3 

2 8 • 3•11 

2 7 - 3 

2® ■ 3 

(-1,4) 

(-1,11) 

(0,1) 

(1,1) 

(1,2) 

(1,4) 

132 

2673 

2 

3 

18 

132 

2 2 • 3 • 11 

3 5 • 11 

2 

3 

2 ■ 3 2 

2 2 • 3■11 

528a2 

4224j 1 

256c2 

384b 1 

96b4 

528a2 

2 4 • 3 • 11 

2 7 ■ 3 - 11 

2 8 

2 7 ■ 3 

2® - 3 

2 4 • 3■11 

(1,11) 

(2,1) 

(2,3) 

(3,1) 

(4,1) 

(5,1) 

2673 

6 

66 

11 

18 

27 

3 5 • 11 

2 ■ 3 

2-3-11 

11 

2 ■ 3 2 

3 3 

4224j 1 

768hl 

8448tl 

1408d2 

768c2 

384g2 

2 7 • 3 • 11 

2 8 - 3 

2 s - 3 - 11 

2 7 • 11 

2 8 - 3 

2 7 - 3 

(5,2) 

(7,4) 

(7, 6) 

(8,1) 

(13,8) 

(14. 1) 

66 

324 

726 

66 

2376 

198 

2-3-11 

2 2 ■ 3 4 

2 ■ 3 - ll 2 

2-3-11 

2 3 ■ 3 3 • 11 

2 - 3 2 • 11 

2112i2 

24a6 


8448c2 

66a2 

8448f2 

2 6 ■ 3 • 11 

2 3 • 3 

2 s - 3 - 11 

2 8 • 3 - 11 

2-3-11 

2 8 • 3 • 11 

(17,2) 

(19, 1) 

(22,1) 

(25,32) 

(31,8) 

(50, 9) 

594 

363 

486 

85536 

8712 

23958 

2 ■ 3 3 ■ 11 

3 ■ ll 2 

2 ■ 3® 

2® - 3® ■ 11 

2 3 - 3 2 ■ ll 2 

2 • 3 2 - ll 3 

2112b2 

4224s2 

768d4 

66c2 

2112r4 

8448q2 

2 6 - 3 - 11 

2 7 - 3 - 11 

2 8 ■ 3 

2-3-11 

2® - 3 - 11 

2 8 • 3 • 11 

(59,16) 

(140,1) 

(155,4) 

(241, 22) 

(695, 8) 


63888 

19602 

96228 

1299078 

3865224 


to 

CO 

CO 

2 ■ 3 4 • ll 2 

2 2 - 3 7 - 11 

2 • 3 10 ■ 11 

2 3 • 3 - 11® 

empty 

2112e4 

8448e2 

2112j2 

2112n4 

2112d4 


2 6 -3 - 11 

2 8 - 3 - 11 

2® - 3 - 11 

2® • 3 • 11 

2® • 3 - 11 
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Table 7: (a, b , c, d) = (0,1,0,1), 6 = -1 • 2 2 , S = {2,3,7,11} 


(-1,1) 

(0,1) 

(1,1) 

2 

1 

2 

2 

1 

2 

128cl 

256c2 

128cl 

2 7 

2 8 

2 7 


Table 8: (a, 6, c, d) = (0,1,0,1), S = -1 • 2 2 , 5 = {2,5,13} 


(-239, 1) 

(-199,32) 

(-83, 64) 

(-63,16) 

(-57,1) 

(-37,16) 

57122 

1300000 

703040 

67600 

3250 

26000 

2 ■ 13 4 

2 s • 5 5 ■ 13 

2 6 ■ 5 ■ 13 3 

2 4 • 5 2 ■ 13 2 

2 • 5 3 - 13 

2 4 • 5 3 ■ 13 

1664k2 

4160i2 

130a4 

4160a4 

8320d2 

4160n2 

2 7 • 13 

2 6 • 5•13 

2-5-13 

2 6 • 5•13 

2 7 ■ 5 - 13 

2 6 • 5 • 13 

(-29,2) 

(-18,1) 

(-12,5) 

(-11,2) 

(-9,13) 

(-8,1) 

1690 

325 

845 

250 

3250 

65 

2 ■ 5 • 13 2 

5 2 ■ 13 

5 - 13 2 

2 ■ 5 3 

2 • 5 3 - 13 

5 • 13 

4160p2 

2080b2 

2080f4 

320c4 

8320al 

2080c2 

2 6 ■ 5 ■ 13 

2 5 • 5•13 

2 5 • 5 • 13 

2 6 -5 

2 7 • 5 - 13 

2 5 • 5•13 

(-7,1) 

(-7,4) 

(-5,1) 

(-3,1) 

(-3,2) 

(-3,4) 


260 

26 

10 

26 

100 

2 • 5 2 

2 2 • 5•13 

2 ■ 13 

2 - 5 

2 ■ 13 

2 2 ■ 5 2 


4160c2 

1664d2 

640a2 

52al 

40a4 

2 7 • 5 

2 6 • 5•13 

2 7 ■ 13 

2 7 ■ 5 

2 2 ■ 13 

2 3 ■ 5 

(-2,1) 

(-1,1) 

(-1,2) 

(-1,5) 

(-1,8) 

(0,1) 

5 

2 

10 

130 

520 

1 

5 

2 

2 - 5 

2-5-13 

2 3 ■ 5 - 13 

1 

160b2 

128cl 

80bl 

8320bl 

65a2 

256c2 

2 5 • 5 

2 7 

2 4 -5 

2 7 ■ 5•13 

5 ■ 13 

2 8 

(1,1) 

(1,2) 

(1,5) 

(1,8) 

(2,1) 

(3,1) 

2 

10 

130 

520 

5 

10 

2 

2 ■ 5 

2-5-13 

2 3 • 5•13 

5 

2 ■ 5 

128cl 

80bl 

8320bl 

65a2 

160b2 

640a2 

2 7 

2 4 -5 

2 7 • 5•13 

5 - 13 

2 5 - 5 

2 7 • 5 

(3,2) 

(3,4) 

(5,1) 

(7,1) 

(7,4) 

(8,1) 

26 

100 

26 


260 

65 

2 ■ 13 

2 2 ■ 5 2 

2 - 13 

2 ■ 5 2 

2 2 ■ 5 - 13 

5 • 13 

52al 

40a4 

1664d2 


4160c2 

2080c2 

2 2 • 13 

2 3 • 5 

2 7 ■ 13 

2 7 ■ 5 

2 6 • 5 - 13 

2 5 • 5 - 13 

(9,13) 

(H,2) 

(12,5) 

(18,1) 

(29, 2) 

(37,16) 

3250 

250 

845 

325 

1690 


2 ■ 5 3 • 13 

2 • 5 3 

5 - 13 2 

5 2 • 13 

2 • 5 • 13 2 

2 4 • 5 3 ■ 13 

8320al 

320c4 

2080f4 

2080b2 

4160p2 


2 7 ■ 5 • 13 

2 6 -5 

2 s - 5 • 13 

2 5 • 5•13 

2 6 ■ 5 - 13 

2 6 • 5 - 13 

(57,1) 

(63,16) 

(83, 64) 

(199,32) 



3250 

67600 

703040 

1300000 



2 ■ 5 3 • 13 

2 4 • 5 2 • 13 2 

2 6 ■ 5 ■ 13 3 

2 5 ■ 5 s ■ 13 

2 ■ 13 4 

empty 

8320d2 

4160a4 

130a4 

4160i2 

1664k2 


2 7 ■ 5 • 13 

2 6 • 5•13 

2-5-13 

2 6 • 5•13 

2 7 ■ 13 
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Table 9: (a, b , c, d) = (0,1,0, -2 ), 6 = 2 3 , S = {2,5,13,17} 


(-4855, 3328) 

(-239,169) 

(-71,8) 

(-37, 26) 

(-33, 20) 

(-31,25) 

4725284096 

-169 

39304 

442 

5780 

-7225 

2 8 • 13 • 17 s 

-1 - 13 2 

2 3 - 17 3 

2-13-17 

2 2 • 5 • 17 2 

-1 - 5 2 • 17 2 

14144v2 

1664h2 

1088c4 

14144s2 

5440h3 

10880f2 

2 6 • 13 • 17 

2 7 ■ 13 

2 6 ■ 17 

2® - 13 ■ 17 

2® - 5 - 17 

2 7 • 5•17 

(-24, 17) 

(-23,16) 

(-7,4) 


(-7,13) 

(-6,1) 

-34 

272 

68 


-3757 

34 

-1-2-17 

2 4 ■ 17 

2 2 ■ 17 


-1 - 13 ■ 17 2 

2 ■ 17 

4352f2 

1088i2 

1088g2 


28288f2 

4352d2 

2 s • 17 

2 6 ■ 17 

2® ■ 17 

2 7 • 5 

2 7 ■ 13 ■ 17 

2 s - 17 

(-5,2) 

(-4,5) 

(-3,2) 

(-2,1) 

(-1.1) 

(0,1) 

34 

-170 

2 

2 

-1 

-2 

2 ■ 17 

-1-2-5-17 

2 

2 

-1 

-1 - 2 

1088j2 

21760q2 

64a2 

256a2 

128c2 

256c2 

2 6 ■ 17 

2 8 - 5 - 17 

2 6 

2 8 

2 7 

2 8 


(2,1) 

(3,2) 

(4,5) 

(5,2) 

(6,1) 


2 

2 

-170 

34 

34 


2 

2 

-1-2-5-17 

2 - 17 

2 ■ 17 


256a2 

64a2 

21760q2 

1088j2 

4352d2 


2 8 

2 6 

2 s • 5 - 17 

2® ■ 17 

2 8 - 17 

(7,4) 

(7,5) 

(7,13) 

(23,16) 

(24,17) 

(31,25) 

68 

-5 

-3757 

272 

-34 

-7225 

2 2 ■ 17 

-1 ■ 5 

-1 ■ 13 ■ 17 2 

2 4 • 17 

-1 • 2 ■ 17 

-1 - 5 2 ■ 17 2 

1088g2 

640d2 

28288f2 


4352f2 

10880f2 

2 6 ■ 17 

2 7 -5 

2 7 • 13 ■ 17 

2® • 17 

2 s ■ 17 

2 7 • 5•17 

(33, 20) 

(37, 26) 

(71,8) 

(239,169) 

(4855, 3328) 


5780 

442 

39304 

-169 

4725284096 


2 2 - 5 • 17 2 

2-13-17 

2 3 ■ 17 3 

-1 ■ 13 2 

2 8 ■ 13 - 17 s 

empty 

5440h3 

14144s2 

1088c4 

1664h2 

14144v2 


2 6 ■ 5 - 17 

2 6 ■ 13 - 17 

2® ■ 17 

2 7 • 13 

2® • 13 • 17 
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Table 10: (a, b, c, d) = (0,1,0, -2), S = 2 3 , S = {2, 7, 29} 


(-181, 128) 
-896 
-1 • 2 7 • 7 
448c6 

2 6 ■ 7 

(-163,116) 
-39788 
-1 - 2 2 ■ 7 3 ■ 29 
12992bc2 

2® ■ 7 ■ 29 

(-45, 29) 
9947 

7 3 ■ 29 
25984i2 

2 7 ■ 7 ■ 29 

(-41,29) 
-29 
-1 - 29 
3712p2 

2 7 • 29 

(-13,16) 
-5488 
-1 ■ 2 4 - 7 3 
448c4 

2® ■ 7 

(-11,8) 

-56 

-1 ■ 2 3 - 7 
448f2 

2® - 7 

(-10,1) 

(-10,7) 

(-9,4) 

(-5,4) 

(-4, 1) 

(-3,1) 

98 

14 

196 

-28 

14 

7 

2 • 7 2 

2 ■ 7 

2 2 - 7 2 

-1 • 2 2 ■ 7 

2 ■ 7 

7 

1792f2 

1792a2 

448a3 

448h2 

1792e2 

896d2 

2 8 ■ 7 

2 8 ■ 7 

2® ■ 7 

2® ■ 7 

2 8 ■ 7 

2 7 - 7 

(-3,2) 

(-2,1) 

(-1.1) 

(-1,2) 


(1,1) 

2 

2 

-1 

-14 


-1 

2 

2 

-1 

-1-2-7 


-1 

64a2 

256a2 

128c2 

448d2 


128c2 

2 6 

2 8 

2 7 

2® ■ 7 

2 8 

2 7 

(1.2) 

(2,1) 

(3,1) 

(3,2) 

(4,1) 

(5,4) 

-14 

2 

7 

2 

14 

-28 

-1-2-7 

2 

7 

2 

2 ■ 7 

-1 ■ 2 2 - 7 

448d2 

256a2 

896d2 

64a2 

1792e2 

448h2 

2 6 ■ 7 

2 8 

2 7 ■ 7 

2® 

2 8 ■ 7 

2® - 7 

(9,4) 

(10, 1) 

(10,7) 

(11,8) 

(13, 16) 

(41, 29) 

196 

98 

14 

-56 

-5488 

-29 

2 2 ■ 7 2 

2 ■ 7 2 

2-7 

-1 • 2 3 ■ 7 

-1 ■ 2 4 - 7 3 

-1 - 29 

448a3 

1792f2 

1792a2 

448f2 

448c4 

3712p2 

2 6 ■ 7 

2 8 ■ 7 

2 s ■ 7 

2®-7 

2® ■ 7 

2 7 • 29 

(45, 29) 

9947 

7 3 ■ 29 

(163,116) 
-39788 
-1 - 2 2 ■ 7 3 - 29 

(181.128) 
-896 
-1 ■ 2 7 - 7 

empty 

empty 

empty 

25984i2 

2 7 - 7 - 29 

12992bc2 

2® - 7 - 29 

448c6 

2® ■ 7 
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Table 11: (a, b, c, d) = (0,1,0, -3), 5 = 2 2 • 3, S = {2,3, 7,11} 


(-122, 9) 

(-111, 64) 

(-97, 56) 

(-85, 49) 

(-69, 16) 

(-53, 7) 

131769 

2112 

56 

1078 

63888 

18634 

3 2 • ll 4 

2 6 • 3 • 11 

2 3 • 7 

2 • 7 2 • 11 

2 4 ■ 3 • ll 3 

2 • 7 • ll 3 

3168z3 

6336b4 

4032bm5 

88704z2 

6336bi4 

88704cs2 

2 5 • 3 2 • 11 

2 6 • 3 2 • 11 

2 6 • 3 2 • 7 

2 7 • 3 2 • 7 • 11 

2 6 • 3 2 • 11 

2 7 • 3 2 - 7 • 11 

(-47, 27) 

(-43, 24) 

(-38, 21) 

(-31, 18) 

(-27, 1) 

(-19, 11) 

594 

2904 

2541 

-198 

726 

-22 

2 • 3 3 • 11 

2 3 • 3 • ll 2 

3 • 7 • ll 2 

-1 • 2 • 3 2 • 11 

2 • 3 • ll 2 

-1-2-11 

12672i2 

6336bd3 

22176f4 

6336m2 

12672c2 

12672bc2 

2 7 • 3 2 • 11 

2 6 • 3 2 • 11 

2 5 • 3 2 ■ 7 • 11 

2 6 • 3 2 ■ 11 

2 7 • 3 2 • 11 

2 7 • 3 2 • 11 

(-15, 8) 

(-15, 14) 

(-13,4) 

(-13,7) 

(-12, 7) 

^ (-9, 4) 

264 

-5082 

484 

154 

-21 

132 

2 3 • 3 • 11 

-1 • 2 • 3 • 7• ll 2 

2 2 • ll 2 

2-7-11 

-1-3-7 

2 2 • 3 • 11 

6336e2 

44352c2 

6336ck3 

88704be2 

2016b2 

6336bh2 

2 6 • 3 2 • 11 

2 6 • 3 2 • 7 • 11 

2 6 • 3 2 - 11 

2 7 • 3 2 • 7 • 11 

2 5 • 3 2 • 7 

2 6 • 3 2 • 11 

(-9, 7) 

( — 7, 3) 

(-7,4) 

(-6, 1) 

(-5, 1) 

(-5,3) 

-462 

66 

4 

33 

22 

-6 

-1 • 2 • 3 • 7 • 11 

2-3-11 


3 ■ 11 

2 • 11 

-1-2-3 

88704bk2 

12672m2 


3168b2 

12672k2 

1152k2 

2 7 • 3 2 • 7 • 11 

2 7 • 3 2 • 11 


2 5 ■ 3 2 ■ 11 

2 7 • 3 2 • 11 

2 7 • 3 2 

(-4.3) 

r (-3, i) 

(-3,2) 

(-2, 1) 

(-1,1) 

(-1,2) 

-33 

6 

-6 

1 

-2 

-22 

-1-3-11 

2 • 3 

-1-2-3 

1 

-1 • 2 

-1-2-11 

3168v2 

1152o2 

576a4 

288c2 

1152r2 

6336ca2 

2 5 • 3 2 • 11 

2 7 • 3 2 

2 6 • 3 2 

2 5 • 3 2 

2 7 • 3 2 

2 6 • 3 2 • 11 

(-1,9) 

(0, 1) 

(1,1) 

(1,2) 

(1,9) 


-2178 

-3 

-2 

-22 

-2178 


-1 • 2 • 3 2 ■ ll 2 

-1 • 3 

-1 ■ 2 

-1-2-11 

-1 • 2 • 3 2 • ll 2 


12672x2 

2304j2 

1152r2 

6336ca2 

12672x2 


2 7 ■ 3 2 • 11 

2 8 • 3 2 

2 7 • 3 2 

2 6 • 3 2 • 11 

2 7 • 3 2 • 11 


(3, 1) 

(3,2) 

(4,3) 

(5, 1) 

(5,3) 

(6, 1) 

6 

-6 

-33 

22 

-6 

33 

2 • 3 

-1-2-3 

-1-3-11 

2 ■ 11 

-1-2-3 

3 • 11 

1152o2 

576a4 

3168v2 

12672k2 

1152k2 

3168b2 

2 7 • 3 2 

2 6 • 3 2 

2 5 • 3 2 • 11 

2 7 • 3 2 • 11 

2 7 • 3 2 

2 5 • 3 2 • 11 

(7,3) 



(9,7) 

(12, 7) 

[13741 

66 



-462 

-21 

484 

2-3-11 



-1 • 2 • 3 • 7 • 11 

-1-3-7 

2 2 • ll 2 

12672m2 



88704bk2 

2016b2 

6336ck3 

2 7 ■ 3 2 • 11 



2 7 • 3 2 • 7 • 11 

2 5 • 3 2 • 7 

2 6 • 3 2 • 11 

(13,7) 

(15,8) 

(15, 14) 

(19,11) 

[27TT1 

(31, 18) 

154 

264 

-5082 

-22 

726 

-198 

2-7-11 

2 3 • 3 • 11 

-1 • 2 • 3 • 7 • ll 2 

-1-2-11 

2 • 3 • ll 2 

-1 • 2 • 3 2 • 11 

88704be2 

6336e2 

44352c2 

12672bc2 

12672c2 

6336m2 

2 7 • 3 2 • 7 • 11 

2 6 • 3 2 • 11 

2 6 • 3 2 ■ 7 • 11 

2 7 ■ 3 2 • 11 

2 7 • 3 2 • 11 

2 6 • 3 2 • 11 

(38, 21) 

(43, 24) 

(47, 27) 

(53, 7) 

(69, 16) 

(85, 49) 

2541 

2904 

594 

18634 

63888 

1078 

3 • 7 ■ ll 2 

2 3 • 3 • ll 2 

2 ■ 3 3 • 11 

2 • 7 ■ ll 3 

2 4 • 3 • ll 3 

2 • 7 2 • 11 

22176f4 

6336bd3 

12672i2 

88704cs2 

6336bi4 

88704z2 

2 5 • 3 2 • 7 • 11 

2 6 • 3 2 • 11 

2 7 • 3 2 • 11 

2 7 • 3 2 • 7 • 11 

2 6 • 3 2 • 11 

2 7 • 3 2 • 7 • 11 

(97, 56) 

(111, 64) 

(122, 9) 




56 

2112 

131769 




2 3 • 7 

2 6 • 3 • 11 

3 2 • ll 4 



empty 

4032bm5 

6336b4 

3168z3 




2 6 • 3 2 • 7 

2 6 • 3 2 • 11 

2 5 • 3 2 • 11 





Table 12: (a, b, c, d) = (0,1,0, -7), 5 = 2 2 • 7, S = {2,5, 7,11} 


(-21,8) 

(-3,1) 

(0,1) 

(3,1) 

(21,8) 

-56 

2 

-7 

2 

-56 

-1 • 2 3 • 7 

2 

-1 • 7 

2 

— 1 • 2 3 • 7 

3136r4 

6272b2 

12544a2 

6272b2 

3136r4 

2 6 • 7 2 

2 7 • 7 2 

to 

00 

“<! 

bO 

2 7 • 7 2 

2 6 • 7 2 
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Table 13: (a, b, c, d) = (1, -1, -4, -1), 5 = 13 2 , S = {2,5,13} 


(-157, 114) 

(-43,157) 

(-17,13) 

(-17,14) 

(-11,8) 

(-9,7) 

65 

-65 

625 

1625 

5 

125 

5 ■ 13 

-1 ■ 5 • 13 

5 4 

5 3 • 13 

5 

5 3 

54080bk2 

54080bk2 

54080dil 

54080bkl 

54080czl 

54080bl2 

2 6 • 5 • 13 2 

2® • 5 • 13 2 

2® • 5 ■ 13 2 

2® ■ 5 ■ 13 2 

2® • 5 ■ 13 2 

2® • 5 • 13 2 

(-7,5) 

(-4,1) 

(-4,3) 

(-4,17) 

(-3,1) 

(-3,2) 

-13 

-65 

5 

-625 

-25 

-5 

-1 ■ 13 

-1 • 5 • 13 

5 

— 1 • 5 4 

-1 ■ 5 2 

-1 ■ 5 

10816bil 

54080cdl 

54080cbl 

54080dil 

54080dcl 

54080cxl 

2® • 13 2 

2® ■ 5 ■ 13 2 

2® • 5 ■ 13 2 

2® ■ 5 ■ 13 2 

2® • 5 ■ 13 2 

2® ■ 5 • 13 2 

(-3,4) 

(-3,11) 

(-3,17) 

(-2,1) 

(-2,3) 

(-2,7) 

65 

-5 

-1625 

-5 

25 

13 

5 ■ 13 

-1 • 5 

— 1 • 5 3 ■ 13 

-1 ■ 5 

5 2 

13 

54080c.dl 

54080czl 

54080bkl 

54080col 

54080dcl 


2® • 5 • 13 2 

2® • 5 • 13 2 

2® • 5 ■ 13 2 

2® • 5 ■ 13 2 

2® ■ 5 ■ 13 2 

2® ■ 13 2 

(-2,9) 

(-1,1) 

(-1,2) 

(-1,3) 

(-1,4) 

(0,1) 

-125 

1 

5 

5 

-5 

-1 

-1 • 5 3 

1 

5 

5 

-1 • 5 

-1 

54080bl2 

10816bel 

54080col 

54080cxl 

54080cbl 

10816bel 

2® • 5 • 13 2 

2® ■ 13 2 

2® • 5 ■ 13 2 

2® ■ 5 ■ 13 2 

2® • 5 ■ 13 2 

2® ■ 13 2 

(1,1) 

(1,2) 

(1,3) 

(2,1) 

(3,1) 

(5,2) 

-5 

-25 

-65 

-5 

5 

-13 

-1 ■ 5 

-1 ■ 5 2 

-1 ■ 5 • 13 

-1 • 5 

5 

-1 • 13 

54080c.ol 

54080dcl 

54080cdl 

54080cxl 

54080cbl 

10816bil 

2® • 5 • 13 2 

2® • 5 • 13 2 

2® • 5 ■ 13 2 

2® • 5 • 13 2 

2® • 5 ■ 13 2 

2® • 13 2 

(7,2) 

(8,3) 

(13,4) 

(14,3) 

(114,43) 


125 

5 

625 

1625 

65 


5 3 

5 

5 4 

5 3 • 13 

5 ■ 13 

empty 

54080bl2 

54080czl 

54080dil 

54080bkl 

54080bk2 


2® • 5 • 13 2 

2® ■ 5 • 13 2 

2® • 5 ■ 13 2 

2® ■ 5 ■ 13 2 

2® • 5 ■ 13 2 
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Table 14: (a, 6, c, d) = (1, -1, -2, -2), S = -1 • 2 3 • 19, S = {2,5,19} 


(-81, 13) 
-593750 
-1 ■ 2 • 5 6 • 19 
115520ckl 

2 6 • 5 • 19 2 

(-43, 89) 
-972800 
-1 • 2 11 • 5 2 ■ 19 
3610il 

2 - 5 • 19 2 

(-17, 16) 

-9025 

-1 • 5 2 • 19 2 
12160il 

2 7 • 5•19 

(-8,9) 

-1250 
-1 • 2 • 5 4 
24320sl 

2 8 • 5 • 19 

-380 

-1 • 2 2 - 5 • 19 
115520ql 

2 6 • 5 • 19 2 

(-3, 1) 

-32 
-1 • 2 5 

23104rl 

2 6 • 19 2 

(-3,4) 

(-2, 1) 

(-2,3) 

(-1. 1) 

(-1,3) 

(-1.H) 

-95 

-10 

-38 

-2 

-40 

-2432 

-1-5-19 

-1-2-5 

-1-2-19 

-1 ■ 2 

-1 • 2 3 • 5 

05 

M 

1 

231040ul 

24320jl 

4864pl 

608fl 

190bl 

38al 

2 7 • 5 • 19 2 

2 8 • 5 • 19 

2 8 • 19 

2 5 - 19 

2-5-19 

2 • 19 

(0, 1) 

(1.1) 

(1.2) 

(1.5) 

(2, 1) 

(3, 1) 

-2 

-4 

-25 

-304 

-2 

10 

-1 • 2 

-1 • 2 2 

-1 • 5 2 

-1 • 2 4 • 19 

-1 • 2 

2 • 5 

4864j 1 

23104bul 

231040bxl 

23104bkl 

92416wl 

115520bxl 

2 8 • 19 

2 6 • 19 2 

2 7 • 5 • 19 2 

2 6 • 19 2 

2 8 • 19 2 

2 6 • 5 • 19 2 

(4, 1) 

(5,2) 

(7,3) 

(9, 29) 

(11. 7) 

(13, 1) 

38 

19 

16 

-65536 

-1280 

2000 

2 • 19 

19 

2 4 

-1-2 16 

-1 • 2 8 • 5 

2 4 • 5 3 

92416hl 

46208bl 

2310411 

23104bt3 

115520bhl 

115520by1 

2 8 • 19 2 

2 7 • 19 2 

2 6 • 19 2 

2 6 • 19 2 

2 6 • 5 • 19 2 

2 6 • 5 • 19 2 

(13791 

(TeTTl 

(T7791 

(18, 11) 

(25, 11) 

(93,41) 

-2888 

50 

-1900 

-4750 

38 

-760 

-1 • 2 3 ■ 19 2 

2 • 5 2 

-1 • 2 2 • 5 2 • 19 

-1 • 2 • 5 3 • 19 

2 • 19 

-1 • 2 3 • 5 • 19 

23104bs2 

24320gl 

115520bml 

24320pl 

23104bjl 

115520by2 

2 6 • 19 2 

2 8 • 5 • 19 

2 6 • 5 • 19 2 

2 8 • 5 • 19 

2 6 • 19 2 

2 6 • 5 • 19 2 

(376, 177) 
-6516050 
-1 • 2 • 5 2 • 19 4 

empty 

empty 

empty 

empty 

empty 

24320hl 

2 8 • 5 ■ 19 
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Table 15: (a, b, c, d) = (1,0,0,1), 6 = -1 • 3 3 , S = {2,3,5} 


(0,1) 

(1,1) 

(1,2) 

(2,1) 



1 

2 

9 

9 



1 

2 

3 2 

3 2 

empty 

empty 

15552b2 

2304j2 

48a4 

48a4 



2 6 • 3 5 

2 8 • 3 2 

2 4 • 3 

2 4 -3 




Table 16: (a, b, c, d) = (1,0, 0, 2), 5 = -1 • 2 2 • 3 3 , S = {2,3, 5} 


(-37, 29) 
-1875 
— 1 • 3 • 5 4 
4320cl 

2 5 • 3 3 • 5 

(-5,4) 

3 

3 

432a4 

2 4 • 3 3 

(-4,3) 

-10 

-1-2-5 

8640cbl 

2 6 • 3 3 • 5 

(-3,1) 
-25 
-1 • 5 2 
4320jl 

2 5 • 3 3 • 5 

(-2,1) 

-6 

-1-2-3 

172811 

2 6 • 3 3 

(-1,1) 

1 

1 

86411 

2 5 • 3 3 

(-1,2) 

(0,1) 

(1,1) 

(2,1) 

(4,7) 


15 

2 

3 

10 

750 


3-5 

2 

3 

2-5 

2 • 3 • 5 3 

empty 

2160wl 

15552b2 

1728hl 

1080jl 

540e2 


2 4 • 3 3 • 5 

2 6 • 3 5 

2 6 • 3 3 

2 3 • 3 3 • 5 

2 2 • 3 3 • 5 



Table 17: (a, 6, c, d) = (1,0,0, -2), 6 = -1 • 2 2 • 3 3 , S = {2, 3, 5} 


(-4,7) 

-750 

-1-2-3-5 3 

540e2 

2 2 • 3 3 • 5 

(-2,1) 

-10 

-1-2-5 

1080jl 

2 3 • 3 3 • 5 

(-1,1) 

-3 

-1-3 

1728hl 

2 6 • 3 3 

(0,1) 

-2 

-1-2 

15552b2 

2 6 • 3 5 

(1,1) 

-1 

-1 

864il 

2 5 • 3 3 

(1,2) 

-15 

-1-3-5 

2160wl 

2 4 • 3 3 • 5 

(2,1) 

(3,1) 

(4,3) 

(5,4) 

(37,29) 


6 

25 

10 

-3 

1875 


2-3 

5 2 

2-5 

-1-3 

3 • 5 4 

empty 

172811 

4320jl 

8640cbl 

432a4 

4320cl 


2 6 • 3 3 

2 5 • 3 3 • 5 

2 6 • 3 3 • 5 

2 4 • 3 3 

2 5 • 3 3 • 5 
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Table 18: (a, b, c, d) = (1,0,0, -3), 5 = -1 • 3 5 , S = {2,3,5} 


(-21,17) 

(-5,1) 

(-3,1) 

(-1,1) 

(-1,2) 


-24000 

-128 

-30 

-4 

-25 

-3 

-1 • 2 6 • 3 • 5 3 

— 1 • 2 7 

-1-2-3-5 

-1 • 2 2 

— 1 • 5 2 

-1-3 

2430gl 

486al 

38880ol 

1944jl 

19440ml 

15552b2 

2 • 3 5 • 5 

2 • 3 5 

2 5 • 3 5 • 5 

2 3 • 3 5 

2 4 • 3 5 • 5 

2 6 • 3 5 

(1,1) 

(1,3) 

(2,1) 

(3,1) 

(3,2) 

(7,5) 

-2 

-80 

5 

24 

3 

-32 

-1-2 

-1 • 2 4 • 5 

5 

2 3 -3 

3 

-1 -2 5 

15552kl 

77760ql 

77760yl 

15552bwl 

15552bdl 

15552bo2 

2 6 • 3 5 

2 6 • 3 5 • 5 

2 6 • 3 5 • 5 

2 6 • 3 5 

2 6 • 3 5 

2 6 • 3 5 

(9,7) 

(11,3) 

(13,9) 

(33,19) 



-300 

1250 

10 

15360 



— 1 • 2 2 • 3 • 5 2 

2 • 5 4 

2-5 

2 10 • 3 • 5 

empty 

empty 

77760vl 

77760bvl 

77760csl 

77760b2 



2 6 • 3 5 • 5 

2 6 • 3 5 • 5 

2 6 • 3 s • 5 

2 6 • 3 5 • 5 
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